1 Squeezing and feedback 



H.M. Wiseman 

Department of Physics, University of Queensland, Queensland 4072 Austràlia 
and 

School of Science, Grifïith University, Nathan, Queensland 4111 Austràlia ** 



Abstract. Electro-optical feedback has many features in common with optical 
nonlinearities and hence is relevant to the generation of squeezing. First, I discuss 
theoretical and experimental results for traveling-wave feedback, emphasizing how 
the "in-loop" squeezing (also known as "squashing") differs from free squeezing. 
Although such feedback, based on ordinary (demolition) photodetection cannot 
create free squeezing, it can be used to manipulate it. Then I treat feedback based on 
noniinear quantum optical measurements (of which non-demolition measurements 
are one example). These are able to produce free squeezing, as shown in a number 
of experiments. FoUowing that I discuss theories showing that intracavity squeezing 
can be increased using ordinary feedback, and produced using QND-based feedback. 
Finally, I return to "squashed" fields and present recent results showing that the 
reduced in-loop fluctuations can suppress atòmic decay in a manner analogous to 
the effect for squeezed fields. 



1.1 Introduction 

In its broadest conception, feedback could be defined to be any mechanism by 
which a system acts upon itself, via an intermediate system. This definition 
would classify, for example, a noniinear refractive index as feedback on a beam 
of light. The light polarizes the medium in which it is propagating, which then 
affects the propagation of the light in that medium. If one is interested in the 
light as a quantum system, then such feedback can be modeled by modifying 
the Hamiltonian for the field. This is a well-known mechanism for gcncrating 
squeezed states of light P . 

In this chapter I am concerned with a different concept of feedback, in 
which the intermediate system is external to the system of interest. That 
is to say, the system is an open system, in constant interaction with its en- 
vironment. Feedback will occur if the change in the environment due to the 
system is significant in affecting the system's dynamics. Since an environment 
is by definition large compared to the system, it is usual for this feedback 
to be ignored. This is the essence of the Born, or perturbative, approach 
to open Systems j^- However, the system's environment can be deliberately 
engineered so that the feedback is important. 
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One obvious way to achieve this is for the environment to include a mea- 
surement apparatus which detects the influence of the system on its sur- 
roundings, a device to aniphfy this measurement, and a mechanism by which 
this amphfied signal controls the dynamics of the system. It is also possible 
to engineer a feedback mechanism which does not involve a measurement 
device, but rather some more direct form of back-coupling from the environ- 
ment to the system. In classical mechanics it is a moot point whether a device 
is designated a "measurement apparatus" . However in quantum mechanics 
the special role of measurement implies that the direct back-coupling may be 
quite distinct from feedback via measurement This is one reason why a 
peculiarly quantum theory of feedback is necessary, and interesting. In this 
chapter I will be concerned only with measurement-based feedback. 

The history of feedback in quantum òptics goes back to the observation of 
sub-shot-noise fiuctuations in an in-loop photocurrent in the mid 1980's by 
two groups |4l5j . The theory of this phenomenon was soon addressed by Ya- 
mamoto and co-workers |til7| . and Shapiro eí aZ The central question they 
were addressing was whether this feedback was producing real squeezing, a 
question whose answer is not as straightforward as might be thought. These 
treatments were based in the Heisenberg picture and used quantum Langevin 
equations where necessary to describe the evolution of system operators. 
They treated the quantum noise only within a linearized approximation. Al- 
though this approximation is probably vàlid for all quantum optical feedback 
experiments performed so far, it would not be vàlid in the "deep quantum" 
regime involving few photons and non-perturbative couplings, as is being 
explored in the so-called "cavity quantum electrodynamics" experiments [Hj . 

More recently an alternative approach to quantum feedback has been 
proposed by myself and Milburn llOIllj , and developed fully by myself ^21 ■ 
This is based on the theory of quantum trajectòries [13114115] . which is an 
application of quantum measurement theory to continuously monitored open 
quantum systems. By treating the measurement explicitly, this theory trans- 
lates the quantum noise of the bath into classical noise in the record of detec- 
tions. It can be shown to be equivalent to an exact (unlinearized) quantum 
Langevin treatment |12j . The advantage of the quantum trajectory method is 
that it allows arbitrary feedback to be treated by the theory, at least numer- 
ically. A particular limit of interest is that of Markovian feedback, in which 
the feedback dynamics can be modeled using a master equation. This result 
was not obtained by the authors using the quantum Langevin treatment. 

A third approach to feedback in quantum òptics is to use the Glauber- 
Sudarshan P function |17ll8ll9j . a quasi-probability distribution. In this the- 
ory, the fields are given an essentially classical description, but negative prob- 
abilities are allowed in order to take into account quantum correlations 
This theory is just as easy to use as the quantum Langevin or quantum 
trajectory theories when the system dynamics which can be linearized. How- 
ever, like the quantum Langevin approach, it is usually intractable when the 
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linearization approximation cannot be made. I will not discuss this theory 
further in this chapter. 

A fourth approach is to treat the electromagnètic field as a stream of 
point-like partides (photons) traveling at the speed of hght. This is essen- 
tially a classical approach, which cannot describe phase properties of the 
fields, but which is adequate if one is interested only in intensity statistics. 
Formally, the in-loop photon arrivals become a self-excited classical point pro- 
cess. This theory was used by Shapiro et al. 8 in addition to their quantum 
operator theory. Similar ideas have subsequently been used by other authors 
[ 2012 l·| . Like the other three approaches mentioned above this approach is 
easily applicable to linearized systems, but unlike them it does not give a full 
quantum description of the in-loop field. Again, I will not discuss this theory 
further in this chapter. 

In the remainder of this chapter I have alternated 'theory' sections, which 
introduce the mathematical apparatus necessary for describing quantum feed- 
back, with 'application' sections, which use the theory to investigate squeez- 
ing, and, where appropriate, discuss experimental results. In both of these 
parallel streams the material is presented in roughly the order in which it 
was developed, but the two streams are not synchronous. 

First I introduce continuum fields, and then show how linearization allows 
feedback onto those fields to be treated analytically, yielding noise spectra 
for in-loop and out-of-loop measurements. Next I discuss the interaction of 
continuum fields with a localized quantum system, giving rise to quantum 
Langevin equations for system operators. This theory is used to describe 
nonlinear measurements (such as QND measurements) of continuum fields, 
and feedback based on the results of these measurements. An alternative to 
the quantum Langevin description is one based on quantum trajectòries. This 
is most useful for illuminating feedback onto the localized systems, and I use 
it to investigate intracavity squeezing. In the Markovian limit the quantum 
trajectory picture of feedback allows one to derive a feedback master equation. 
This is of most use for describing dynamics which cannot be linearized, such 
as that of a strongly driven two-level atom. This turns out to be precisely 
what is needed to revisit the question of in-loop squeezing in terms of what 
the atom 'sees'. 

1.2 Continuum Fields 
1.2.1 Canonical Quantization 

Let the fundamental field be the vector potential A(r,t) in the Coulomb 
gauge 



V • A(r,í) = 0. 
The free Lagrangian density for this field is [23 



(1.1) 




(1.2) 
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where Eq is the permittivity of free space and c is the speed of light. The 
elèctric E and magnètic B fields are defined by 

E = -À ; B = V X A. (1.3) 

From H1.2|l . the canonical field to A is — eoE. Thus, in quantizing the field, 
these obey the canonical commutation relations 

[Mr,t)^Ekir',t)] = -^^S,,Siir-r'). (1.4) 

Here denotes a three-dimensional transverse Dirac delta-function, which 
is necessary to be compatible with the constraint of (|l.l|l Z^ . Note that 
the Heisenberg picture operators in the canonical commutation relations are 
at equal times. In the Schròdinger picture, the same relations hold, but the 
time argument is omitted. The Euler-Lagrange (which is also the Heisenberg) 
equation of motion from (|1.2|l is the wave equation 

À = c^V^A. (1.5) 

Now consider the case of a beam of polarized light. That is to say, consider 
only one component A of A and let its spatial variation be confined to one 
direction, say z. This simplifies the analysis, and is also appropriate for de- 
termining the inputs and outputs of a quantum optical cavity. In reality, the 
transverse spatial extent of the beam would be confined to some area A which 
is determined by the area of the optical components involved However, 
as long as the x and y extensions are much greater than a wavelength, the 
beam can be approximated by plane waves. The appropriate wave equation 
is 

À = c^a^A, (1.6) 

of which I am interested only in the forward propagating solutions 

A{z,t + t')=A{z-ct',t). (1.7) 

If the field is reflected off a cavity mirror (say at z = 0) thcn the direction 
of z will change at the point of reflection. This is why only one direction of 
propagation need be considered. The field for z < is incoming and that for 
> is outgoing. The canonical commutation relation is now 

[A{z,t),E{z',t)] = -t^S{z-z'). (1.8) 

EqA 

Solutions for A and E satisfying the wave equation H1.6|l can be con- 
structed using the annihilation and creation operators for the modes of fre- 
quency u, which satisfy 

[aitü),aHLü')] = S{Lü~Lü'). (1.9) 
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They are 



E{z) 



A{z) 




This expression for the fields in terms of annihilation and creation opera- 
tors for a continuum of modes defines the sense in which thcy are composed 
of photons of definite frequency. However, this sense is quite unhke the naive 
picture of a beam of hght made up of (possibly differcnt freqüències of) 
photons, hm'thng through space at the speed of hght. Each mode is spread 
over all space, so there is no way in which a photon, as an excitation of 
such a mode, can move at all. To define an annihilation operator 6(2, t) for 
a localized photon of a particular frequency, it would be necessary to sum 
many different mode operators. Such operators can be defined, with slight 
variations in the detalls of the definition |14I23| . The various definitions are 
effectively equivalent in application to quantum optical problems. The au- 
thors of (141231 construct the localized annihilation operator from the mode 
annihilation operators 0(0;). Here, just for variation, I am introducing a dif- 
ferent definition for b{z, í), constructcd from the original fields in space-time, 
A{z,t) and E{z,t). 

As established above, A and —EqE are canonically conjugate variables at 
each point in space-time. Motivated by the analogy with position and mo- 
mentum, a local annihilation operator for an oscillator of angular frequency 
ivo can be defined as 



If the beam contains only photons of a frequency near ujq then it is apparent 
from itTTÏÏll that we can approximate b{z,t) by 




E{z,t) .(1.12) 



In terms of the mode operators, b{z,t) is given by 




(1.13) 




(1.14) 



From this we can calculate 




(1.15) 



cS{z - z') 



(1.16) 
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which explains the choice of normalization in (|1.12(l . Note that the approxi- 
mations involved here apply only if all of the light is at a frequency close to 
üjQ. Thus, the width of the Dirac S function in H1.16|l should be understood 
to be much greater than a wavelength of light 27rc/a;o· 

The rotating exponential exp[ia;o(í — -^/c)] means that a "localized pho- 
ton" of frequency lüq has a slowly varying annihilation operator b{z,t). This 
operator also has the same property as the vector potential (|1.7(l . obeying 

biz,t + t') ^b{z- ct',t) (1.17) 

in free space. If only freqüències near ujq are significantly excited, then the 
time-flux of energy can be easily seen to be 

W{z,t) c:ihLüob^z,t)b{z,t). (1.18) 

Thus, the annihilation operator ò(z, t) conforms to one's naive expectations, 
with òt(z, t)b{z, t) being the photon flux (photons per unit time) passing z at 
time t. 



1.2.2 Photodetection 

From the above discussion it should be apparent that it is not sensible to 
talk about a photodctector for photons of frequency which has a response 
time comparable to or smaller than uJq^ ■ Therefore for practical purposes a 
photodctector is equivalent to an energy flux meter. In either long 
as we are not interested in times comparable to uj^^ ^ we can assume that 
the signal produced by an ideal photodctector at position zi is given by the 
operator 

/(í) =òí(í)òi(í), (1-19) 

where bi{t) = b{zi,t). Here I have ignored any factors of elèctric charge etc. 
which are sometimes included but which are actually nominal. 

In experiments involving làsers, it is often the case (or at least it is harm- 
less to assume |21]) that bi{t) has a mean amplitude /3 = (òi(í)). Without 
loss of generality, I will take P to be real. In all that foUows I will also assume 
that we are considering stationary statistics. That is, we are taking the long 
time limit of a system with a stationary state. In that case, only if the corre- 
lations of interest in the intensity of the beam of light have a characteristic 
time satisfying 

Tcor » (1.20) 

is it permissible to linearize 11.19|l . This means approximating it by 

I{t)=(3^+6I{t)=l3^+(3Xi{t), (1.21) 
where 

Xi{t) = bi{t) + b\{t) - 2f3 (1.22) 
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is the amplitude quadrature fluctuation operator for the continuum field. For 
the linearization to be vàlid the fluctuations must be small as well as slow: 

+ t)Xi(í)) « for T ~ Teorr. (1.23) 

It is useful also to define the phase quadrature fluctuation operator 
Yi{t) = -tbi{t)+ib\{t). (1.24) 

For free fields, which obey H1.17|l . these obey the commutation relations 
[Xi(í),ri(í')] = 2^<5(í-í'). (1.25) 

If we define the Fourier transformed operator 

/oo 
dtXi{t)e-"^* (1.26) 
-oc- 

and similar ly for ^1(1^) then 

[Xi{lü),Yi{lü')] ^ AniSiuj + uj'). (1.27) 

For stationary statistics as we are considering, {Xi{t)Xi{t')) is a function of 
t — t' only. From this it follows that 

{Xi{lü)Xi{uj')) cx S{lü + Lü'). (1.28) 

Because of the singularities in equations (|1.27|l and H1.28|l . to obtain a 
finite uncertainty relation it is more useful to consider the spectrum 

1 í°° 

Siii^) = 7^ / (Xiiüj)Xi{~u;'))diü' (1.29) 



2n 



00 
iüjt 



e-"^\X{t)X{0))dt = (Xi(w)Xi(O)). (1.30) 

Then it can be shown that for a stationary free field [S], 

Sf{uj)SY{uj) > 1. (1.31) 

From this it is obvious that a coherent continuum field [17118119] is one such 
that for all ui 

S?{lü) = 1, (1.32) 

where Q = X or Y (or any intermediate quadrature). This is known as the 
Standard quantum limit or shot-noise limit. A squeezed continmmr field is one 
such that, for some uj and some Q, 

S^iiü) < 1. (1.33) 
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The physical significance of (uj) is apparent from H1.2l·|l : it can be 
experimentally determined as 

/oo 
*(/(í),/(0))dí, (1.34) 
-oo 

where {A, B) = (AB) — {A){B). In fact it is possible to determine S^{u!) for 
any quadrature Q in a similar way. Putting the field of interest through a low- 
reflectivity beam splitter, while reflecting another field with a large coherent 
amplitude ofF the same beam splitter, a coherent amplitude can be added to 
the beam of interest. If this contribution is sufficiently large it will dominate 
the total coherent amplitude of the beam. Since the added component can 
have any chosen phase with respect to the original beam, the new linearized 
intensity fluctuation operator will be proportional to any chosen quadrature 
fluctuation operator. This technique is known as homodyne detection. In 
practice, balanced homodyne detection using a 50-50 beam splitter and two 
photodetectors is preferable, but the principle is the same j25j . 



1.3 In-loop "Squeezing" 
1.3.1 Description of the Device 

The simplest form of quantum optical feedback is shown in Fig. 11.11 This 
was the scheme considered by Shapiro et al.. In our notation, we begin with 
a field bo = b{zo, t) as shown in the diagram. We will take this field to have 
stationary statistics with mean amplitude (3 and fluctuations 

\[Xo{t)+iYo{t)]=bo{t)-(i (1.35) 

characterized by arbitrary spectra Sq [lo) , Sq {u)) . This field is then passed 
through a beam splitter of transmittance rji (í) . By unitarity, the diminution 
in the transmitted field by a factor y^rji{t) must be accompanied by the 
addition of vacuum noise from the othcr port of the beam splitter T. The 
transmitted field is 

bi{t) = Vmit-T,)bo{t-n) + ^f]iit-n)u{t-Ti). (1.36) 

Here ti — {zi — zo)lc and I am using the notation 

í^EEl-íy. (1.37) 

The operator u{t) represents the vacuum fluctuations. The vacuum is special 
case of a coherent continuum field of vanishing amplitude {v{t)) — 0, and so 
is completely characterized by its spectrum 

S^iu:) = 1. (1.38) 

Since the vacuum fluctuations are uncorrelated with any other field, and have 
stationary statistics, the phase and time arguments for v{t) are arbitrary. 
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Z=Zo Z=Zi 




Fig. 1.1. Diagram for a traveling-wave feedback experiment. Traveling fields are 
denoted 6 and photocurrent fluctuations 51. The first beam splitter transmittance 
rji is variable, the second 772 fixed. The two vacuum field inputs are denoted and 



The beam-splitter transmittance í/i(í) in H1.36|) is time-dependent. This 
time-dependence can be achieved experimentally by a number of means. For 
example, if the incoming beam is eUiptically polarized then an electro-optic 
modulator (a device with a refractive index controUed by a current) will alter 
the orientation of the elhpse. A polarization-sensitive beam sphtter will then 
control the amount of the light which is transmitted, as done, for example, in 
|2t)|. As the reader will no doubt have anticipated, the current used to control 
the electro-optic modulator can be derived from a later detection of the light 
beam, giving rise to feedback. Writing 771 (í) = fïi+ Srji (t) , and assuming that 
the modulation of the transmittance is small {6r]i{t) <C ?7i), one can write 



Vm(t) = Vm + {i/Vm)kSrj,{t). (1.39) 

Continuing to follow the path of the beam in Fig. Il.l·l it now enters a 
second beam-splitter of constant transmittance 772- The transmitted beam 
annihilation operator 



hit) = y/rj2bi(t - T2) + VmtJ.{t - T2), 



(1.40) 
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where T2 — (^2 ~ ^i)/c and /i(í) represents vacuum fluctuations like vit). The 
reflected beam operator is 

bsit) =V^b,{t-T2)- ^H{t-T2). (1.41) 

Using the approximation l|1.39(l . the linearized quadratura fluctuation oper- 
ators for Ò2 are 

X2{t) = VmmMt - T2) + ^/m/mí3Sm{t ~ T2) 

+ Vm^X4t-T2) + V^X^{t-T2), (1.42) 



Y2{t) - v^ro(í - T2) + v^r,(í - T2) + VïhY^it - T2), (1.43) 

where T2 = T2 + ri . Similarly for Ò3 we have 



x-iit) = VmmXoit - T2) + VfkJvi fiSrnit - T2) 

+ V^X,{t-T2)- ^Xf,{t-T2), (1.44) 



Y3{t) = V^Yo{t-T2) + V^Y,{t-T2)-^Y^it-T2). (1.45) 



The mean fields for Ò2 and Ò3 are yj'r\\r\2 P and -/ryi?72 /3 respectively. Thus, 
if these fields are incident upon photodetectors, the respective hnearized pho- 
tocurrent fluctuations are, as explained in Sec. 



5l2{t) = Vm^/3^2(í), (1.46) 
Sh{t) = Vm^í^Xsit). (1.47) 

Here I have assumed perfect efhciency detectors. To model inefhcient detec- 
tors it is necessary to add further beam sphtters, with transmittance equal to 
the detection efficiency, in front of the detectors. The effect of this has been 
considered in detail in |26l27j . 

Having obtained an expression for Sl2{t) we are now in a position to 
follow the next stage in Fig. Il.ll and complete the feedback loop. We set the 
modulation in the transmittance of the first beam-sphtter to be 



00 



'5^1 W = / Ht^hit - TO - t')dt\ (1.48) 

V2P^ Jo 

where g is a dimensionless parameter representing the low-frequency gain of 
the feedback loop. The response of the feedback loop, including the electro- 
optic elements, is assumed to be linear for small fluctuations and is charac- 
terized by the electrònic delay time tq and the response function h{t'), which 
satisfies h{t) = for t < 0, h{t) > for í > and h{t')dt' = 1. 

1.3.2 Stability 

Clearly the feedback can only affect the amplitude quadrature X. Putting 
lfriH|l into itLl^ yiclds 



X2{t) = VmmXoit -T2) + g h{t')X2{t - r - t')dt' 

JO 

+ Vm^X,{t^T2) + V^X^,{t-T2), (1.49) 



1 Squeezing and feedback 11 



where T ^ tq + T2 = tq + ri + T2. This is easy to solve in Fourier space, 
providing that X2 is a stationary stochastic process. This will only be the 
case if the feedback is stable. Using Standard feedback and control theory 
I2H1, the Nyquist stability criterion is 

Re[s] < 0, (1.50) 

where s is any solution of the characteristic equation 

1 - gií(s) exp(-sT) = 0, (1.51) 

where H{s) denotes the Laplace transform H{s) = dte~^'^h{t). 

First I show that a sufficient condition for stabihty is |f;| < 1. Looking for 
instabihty, assume that Re[s] > 0. Thcn 



\H{s) 



'ST\ 



< / dth{t) = 1. (1.52) 



Thus undcr this assumption the characteristic equation cannot be satisfied 
for l^l < 1, so this regime wiU always be stable. If g > 1 then it is not 
difhcult to show that there is a positive s which will solve (|1.51|) . Thus it is 
a necessary condition to have g < 1. lí g < —1, the stability of the feedback 
depends on T and the shape of h{t). However, it turns out that it is possible 
to have arbitrarily large negative low-frequency feedback (that is, g —00), 
for any feedback loop delay T, provided that h{t) is broad enough. The price 
to be paid for strong low-frequency negative feedback is a reduction in the 
bandwidth of the feedback, the width of |/i(ijj)p. 

To see this, consider the simplest smoothing function h{t) — 76^'^*. The 
condition for marginal stability is that there is a solution to (|1.51() for s = iu. 
That is, 

1 = gexpí-iujT)— — . (1.53) 

7 + icj 

For the imaginary part of the right-hand side to vanish wc rcquire 

tancjT = -cj/7. (1.54) 

As we will see, for large \g\ we will require 7 <C T^^ in which case the 
Solutions on H1.54|l can be approximated by a;„ = {2n + \)n/2T. Undcr the 
same approximation wc can ignore 7 compared to lo in lfT35|l to get 

l^~\g\{~^r^'-^^. (1.55) 

Clearly for n odd this cannot be satisfied and so the system will be stable. 
However for n even we have 



'^'(2n+l)^' 



(1.56) 
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which can be satisfied (indicating marginal stability). In order to avoid this 
for all n we require 

where here we see that 7 ^ T^^ for large negative g. Now the bandwidth of 
the feedback is i? ~ 27. Thus we have finally the approximate inequality 

which shows how a finite delay time T and large negative feedback — g S> 1 
reduces the possible bandwidth of the feedback. 



1.3.3 In-loop and Out-of-loop Spectra 

Assuming then that the feedback is stable we can solve H1.49|l for X2 in the 
Fourier domain: 

X2(w) = exp(-zwr2) ~— . . (1.59) 

1 — gn(üj) exp(— ^ti;i j 

From this the amplitude quadrature spectrum is easily found from (|1.28(l and 
to be 

cX/ N VlV2S^ (UJ) + 7]2mS^ {UJ) + fÏ2S^ {cj) 
^2 [^) - 



|1 -5/i(w)exp(-íiwT)|2 
mV2[S^{u;)-l] 



|1 -5/i(w)exp(-iwT)|2 



(1.60) 



From these formulae the effect of feedback is obvious: it multiplies the 
amplitude quadrature spectrum at a given frequency by the factor |1 — 
ghiüj) exp(— ïa;T)|~2. At low freqüències, this factor is simply {\ — g)~^, which 
is why the feedback was classified on this basis into positive {g > 0) and nega- 
tive {g < 0) feedback. The former will increase the noise at low frequency and 
the latter will decrease it. However at higher freqüències, and in particular at 
múltiples of tt/T, the sign of the feedback will reverse and g < will result 
in an increase in noise and vice- versa. This is shown clearly in the theoretical 
investigations of Shapiro et ai. All of these results make perfect sense in the 
context of classical light signals, except that in that case we would not worry 
about vacuum noise. This is equivalent to assuming that the original noise is 
far above the shot-noise limit, so that one can replace 1 + rii'ri2[S^ (uj) — 1] 
by r]iT]2S^ (cu). This gives the result expected from classical signal process- 
ing: the signal is attenuated by the beam splitters and either amplified or 
suppressed by the feedback. 
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The most dramàtic efïect is of course for large negative feedback. For 
sufíiciently large —g it is clear that one can make 

S^iüü) < 1 (1.61) 

for some lu. This effect has been observed experimentally many times with 
different systems involving feedback |4I5I7I29I26I27I3U| . Without a feedback 
loop this sub-shot-noise photocurrent would be seen as evidence for squeezing. 
However, there are a number of reasons to be very cautious about applying 
the word squeezing to this phenomenon. Two of these reasons are theoretical, 
and are discussed in the following two sub-sections. The more practical reason 
relates to the out-of-loop beam Ò3, which 1 will now discuss. 

From H1.44(l . the X quadrature of the beam 63 is, in the Fourier domain, 



Xaiuj) = exp(-ïwr2) ^/tJ^XoÍlü) + Vmvï X^{uj) - .Jrj^X^{L·j) 



+ Vfh/mghiu)eM~^^T)X2{Lü). (1.62) 

Here I have substitutcd for Srji in terms of X2 ■ Now using the above expression 
H1.59|l gives 



X3{uj) = exp{-iujT2) 



y/mvi Xq{uj) + ^/rï^X^{uj) 



1 — gh{u}) exp{—ÍLüT) 
rï^[l - gh{uj) exp{-iujT)/ri2]X^{üj) 



1 — gh{uj) exp(— iwT) 
This yields the spectrum 

^ l+^2r/l(5o^-l) 
^3 [UJ) 



(1.63) 



1 — gh(ijj) exp(— ia-'T)p 
2Kc[gh{L·j) exp(-ít^r)] + g^\h{uj)\'^ /ri2 
|1 -5/i(w)exp(-iwT)|2 



(1.64) 



The denominator is identical to that in the in-loop case, as are the first two 
terms in the numerator. But there are additional terms in the numerator 
which indicate that there is extra noise in the out-of-loop signal. 
The expression (11.64(1 can be rewritten as 

\1- gh(üj)exp{-iujT)\'^ 

From this it is apparent that, unless the initial beam is amplitude squeezed 
(that is, unless Sq{uj) < 1 for some lu) the out-of-loop spectrum will always 
be greater than the shot-noise-limit of unity. In other words, it is not possible 
to extract the apparent squeezing in the feedback loop by using a beam split- 
ter. In fact, in the limit of large negative feedback (which gives the greatest 
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noise reduction in the in-loop signal), the out-of-loop amplitude spectrum 
approaches a constant. Considering a frequency lo such that h{Lü) exp(— icjT) 
is real and positive, one finds that 

lim {(jj) = ri:^^ . (1.66) 

Thus the more hght one attempts to extract from the feedback loop, the 
higher above shot-noise the spectrum becomes. 

This resuh is counter to an intuition based on classical hght signals, where 
the efïect of a beam sphtter is simply to spht a beam so that both outputs 
would have the same statistics. The reason this intuition fails is precisely 
because this is not ah that a beam sphtter does; it also introduces vacuum 
noise which is anticorrelated at the two output ports. The detector for beam 
62 measures the amphtude fluctuations X2, which are a combination of the 
initial fluctuations Xq, and the two vacuum fluctuations X^, and X^. The 
first two of these are commou to the beam 63, but the last, X^, appears 
with opposite sign in X^. As the negative feedback is turned up, the first two 
components are successfuUy suppressed, but the last is actually amplified. 
Note that the result in H1.66|l holds no matter how large So{ío) is compared 
to unity. 



1.3.4 Commutation Relations 

Under normal circumstanccs (without a feedback loop) one would expect a 
sub-shot noise amplitude spectrum to imply a super-shot-noise phase spec- 
trum. However that is not what is found from the theory presented here. 
Rather, the in-loop phase quadrature spectrum is unaffected by the feed- 
back, being equal to 

S^{Lj)^l + mV2[S^{iü)~l]. (1.67) 

It is impossible to measure this spectrum without disturbing the feedback 
loop because all of the light in the 62 beam must be incident upon the pho- 
todetector in order to measure JÍ2. However, it is possible to measure the 
phase-quadrature of the out-of-loop beam by homodyne detection. This was 
done in which verified that this quadrature is also unaffected by the 
feedback, with 

S^icj) = l + mfÍ2K(^)-'^]- (1-68) 

For simplicity, consider the case where the initial beam is coherent with 
S(^{lü) = S^{üj) = 1. Thcn S^{üj) = 1 and 

SI\lü)S^{üj) = |1 - gh{uj) exp{~iujT)\-^. (1.69) 

This can clearly be less than unity. This represents a violation of the uncer- 
tainty relation l|1.31|l which foUows from the commutation relations (ll.27|l . 
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In fact it is easy to show (as done first by Shapiro et al. 8 ) from the solution 
H1.59|l that the commutation relations l|1.27|l are false for the field Ò2 and 
must be replaced by 

1 — gn[Lüj exp(— líjj j 

which explains how H1.69(l is possible. 

At first sight, this apparent violation of the canonical commutation re- 
lations would seem to be a major problem of this theory. In fact, there 
are no violations of the canonical commutation relations. As emphasized in 
Sec. 11.2.11 the canonical commutation relations (|1.8() are between fields at 
different points in space, at the same time. It is only for free fields (traveling 
forward in space for an indefinite time) that one can replace the space differ- 
ence z with a time difference t = z/c. Field Ò3 is such a free field, as it can 
be detected an arbitrarily large distance away from the apparatus. Thus its 
quadratures at a particular point do obey the time-difference commutation 
relations H1.25|l . and the corresponding Fourier domain relations (|1.27() . But 
the field Ò2 cannot travel an indefinite distance before being detected. The 
time from the second beam splitter to the detector T2 is a physical parameter 
in the feedback system. 

For times shorter than the total feedback loop delay T it can be shown 
that 

[X2{t),Y2{t')]^2i5{t-t') ioï\t-t'\<T. (1.71) 

Now the field Ò2 is only in existence for a time T2 before it is detected. Be- 
cause T2 < T, this means that the time-difference commutation relations 
between different parts of field &2 are actually preserved for any time such 
that those parts of the field are in existence, traveling through space towards 
the detector. It is only at times greater than the feedback loop delay time T 
that non-standard commutation relations hold. To summarize, the commu- 
tation relations between any of the fields at different spatial points always 
holds, but there is no reason to expect the time or frequency commutation 
relations to hold for an in-loop field. Without these relations, it is not clear 
how "squeezing" should be defined. Indeed, it has been suggested [31] that 
"squashed light" would be a more appropriate term for in-loop "squeezing" 
because the uncertainty has actually been squashed, rather than squeezed 
out of one quadrature and into another. 



1.3.5 Semiclassical Theory 

A second theorctical reason against the use of the word squeezing to describe 
the sub-shot-noise in-loop amplitude quadrature is that (providing beam 60 
is not squeezed), the entire apparatus can be described semiclassically. In a 
semiclassical description there is no noise except classical noise in the field 
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amplitudes, and shot noise is a result of a quantum detector being driven by a 
classical beam of light. That such a description exists might seem surprising, 
given the importance of vacuum fluctuations in the explanation of the spectra 
in Sec. 11.3.51 However, the semiclassical explanation, as explored in [8i26.32| . 
is even simpler. 

For the fields, let us use the same symbols as before, but with a cl super- 
script to rcmind us that these are classical variables rather than operators. 
Then the only irreducible source of noise in the system is the shot noise at 
the two detectors A: = 2,3. This arises from the assuniption that a classical 
light field of frequency uiq and time-flux of energy induces photo-electron 
emissions as a random process at rate /hujo. Here h appears as a univer- 
sal phenomenological constant relating the output of photodetectors to the 
incoming light. In addition to this irreducible noise, the light field itself may 
have (classical) noise which is represented by the amplitude fluctuation vari- 
able X'^'(í), which is scaled in the same way as the quantum operator X{t) 
has been. 

In the linearized approximation we have been working in, the Poissonian 
shot-noise fluctuations can be approximated by Gaussian fluctuations, giving 



= (1-72) 

where Ik is the mean value of the photocurrent and (t) represents indepen- 
dent white-noise sources obeying 

{^kit%{t))^5,kS{t-t'). (1.73) 



In our case we have ^/T2 = /3-^?yi?72 and vTa = P^/rjï^ 

Using the above expression H1.48|l for Sr], the fluctuation in the field inci- 
dent on detector 2 is 

+ g h{t')[Xf{t-T-t')+Í2{t-T-t')]dt'. (1.74) 



^0 

Assuming stable feedback (the conditions are the same as before) and solving 
this in Fourier domain we get 

^ci^^^) ^ exp(-ii^r2)^/?7^X^'(tj) + gh{uj) g^y>{-íujT%{í^) 

1 — gh{üj) exp(— ia;T) 
But the photocurrent noise itself is proportional to 

X^ (iv) + = j— — — . (1.76) 

1 — gn[Lü) exp(— iciji ) 

From (|1.34(l we find the spectrum 

S^ico)^ l + ry2^i[^o^(a;)-l] 

|l-5/^(w)exp(-^CJ^)|2' 
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where I have defined 



(1.78) 



as the spectrum which would be observed with no feedback and no beam 
splitters. With this Identification, the expression (|1.77|l is identical with H1.6Í)|I 
derived using a quantized field. 

In this formulation, the sub-shot noise of the in-loop current is no surprise 
at aU. It is simply a result of feeding back an amphfied negative version of the 
shot noise randomly produced at the detector so that the current at a later 
time win be anticorrelated with itself. The low noise is in the photocurrent 
only, not in the hght beams which here are all completely classical. This ex- 
plains why the out-of-loop detector does not produce a sub-shot noise signal. 
The photocurrent noise at that detector is proportional to 



Since Ç3 is independent of everything else, it is noise which cannot be reduced 
by the feedback. Only the classical noise Xq' can be reduced, and that at 
the expense of introducing the extra uncorrelated noise ^2- Calculating the 
spectrum Sxi^) again gives the same answer as the quantum theory. 

1.3.6 QND Measurements of In-loop Beams 

On the basis of the above semiclassical theory it might be thought that all of 
the calculations of noise spectra made in this section relate only to the noise 
of photocurrents and say nothing about the noise properties of the light 
beams themselves. While this appears to be the case from the semiclassical 
theory, it is not really true, as can be seen from considering quantum non- 
demolition (QND) measurements 0. In this context, a QND detector is one 
which can measure the intensity of light without absorbing it. Such devices 
cannot be described by semiclassical theory, which shows that this theory is 
not complete and hence cannot be expected to provide the correct intuition 
about the state of the light itself. 

A specific model for a QND device will be considered in Sec. 11.5.11 Here 
we simply assume that a perfect QND device can measure the amplitude 
quadrature X of a continuum field without disturbing it beyond the neces- 
sary back-action from the Heisenberg uncertainty principle. In other words, 
the QND device should give a read-out at time t which can be represented 
by the operator X{t). The correlations of this read-out will thus reproduce 
the correlations of X{t). For a perfect QND measurement of X2 and X3, the 
spectrum will reproduce those of the conventional (demolition) photodetec- 
tors which measure these beams. This confirms that these detectors (assumed 



Ut) + xfit) = ut) + Vm^2X^\t - T2) + 9V^2 



X / h{t')[X:^\t - T - t') + ^2{t - T - t')]dt' . (1.79) 
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perfect) are indeed recording the true quantum fluctuations of the light im- 
pinging upon them. 

What is morè interesting is to consider a QND measurement on Xi . That 
is because the set up in Fig. ll.l·l is equivalent (as mentioned above) to a set up 
without the second beam spUtter, but instead with an in-loop photodetector 
with efficiency 772. In this version, the beams Ò2 and 63 do not physically exist. 
Rather, hi is the in-loop beam and X2 is the operator for the noise in the 
photocurrent produced by the detector. As shown above, this operator can 
have vanishing noise at low freqüències for g —00. However, this is not 
reflected in the noise in the in-loop beam, as recorded by our hypothetical 
QND device. FoUowing the methods of Sec. 11.3.31 the spectrum of Xi is 

gx.. ^ l + Vi[Soiuj)-l]+g^\hiuj)\^fj2/m 
' \l-gh{uj)exp{-iujT)\^ 

In the limit g —>■ —00, this becomes at low freqüències 

^f(0)^i^, (1.81) 

which is not zero for any detection efficiency 772 finitely less than one. Indeed, 
for 7^2 < 0.5 it is above shot-noise. 

The reason that the in-loop amplitude quadrature spectrum is not re- 
duced to zero for large negative feedback is that the feedback loop is feeding 
back noise ^^(í) in the photocurrent fluctuation operator X2(t) which is 
independent of the fluctuations in the amplitude quadrature Xi (í) of the in- 
loop light. The smaller 772 is, the larger the amount of extraneous noise in 
the photocurrent and the larger the noise introduccd into the in-loop light. 
In ordcr to minimize the low-frequency noise in the in-loop light, there is an 
optimal feedback gain. In the case 5o(w) = 1 (a coherent input), this is given 

by 

5opt = --^, (1.82) 

1 - ?72 

giving a minimum in-loop low-frequency noise spectrum 

5f(0)„,in = 1-772. (1.83) 

The fact that the detection efficiency does matter in the attainable in-loop 
squeezing shows that these are true quantum fluctuations. 



1.3.7 A Squeezed Input 

Although the feedback device discussed in this section cannot produce a free 
squeezed beam, it is nevertheless useful for reducing classical noise in the 
output beam 63. It is easy to verify the result of [201 that if one wishes to 
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reduce classical noise Sq {lü) — 1 at a particular frequency íü, then the optimum 
feedback is such that 

gh{Lü)exp{-iujT) = -7nV2[S^ {cj) - 1]. (1.84) 

This gives the lowest noise level in the amphtude of &3 at that frequency 



l + T]2r)i[SQ (UJ) - V 



For large classical noise we have feedback proportional to Sq (uj) and an 
optimal noise value of 1/í72, as this approaches the limit of H1.66|l . The inter- 
esting regime is the opposite one, where Sq (uj) — 1 is small, or even nega- 
tive. This last case corresponds to to a squeezed input beam. Putting squeez- 
ing through a beam splitter reduces the squeezing in both output beams. In 
this case, with no feedback the residual squeezing in beam 63 would be 

53^(w)„o = 1 + mmiSni^) - 1], (1-86) 

which is closer to unity than Sq{uj). The optimal feedback (the purpose of 
which is to reduce noise) is, according to (11.841) . positive. That is to say, 
destabilizing feedback actually puts back in beam Ò3 some of the squeezing 
lost through the beam splitter. Since the required round-loop gain (|1.84|) is 
less than unity, the feedback loop remains stable. 

This result highlights the nonclassical nature of squeezed fluctuations. 
When an amplitude squeezed beam strikes a beam splitter, the intensity at 
one output port is anticorrelated with that at the other, hence the need 
for positive feedback. Preliminary observations of this effect were reported 
in . Of course the feedback can never put more squeezing into the beam 
than was present at the start. That is, S';^(a;)opt always lies between S^ (uj) 
and 5';^(a;)no· However, if we take the limit 771 ^ 1 and Sq (lo) (perfect 
squeezing to begin with) then all of this squeezing can be recovered, for any 
?72 . This effect might even be of practical use for attenuating highly squeezed 
sources, such as those produced by làser diodes, while retaining most of the 
squeezing 



1.4 Quantum Langevin Equations 

In the preceding section, only continuum fields were considered, with all op- 
tical and electro-optical devices being treated as classical (i.e. deterministic) 
Systems. Often one wishes to consider feedback onto quantum systems which 
may introduce extra noise terms into the equations. To do this one needs a 
theory for describing the dynamics of localized quantum systems interact- 
ing with continuum fields. In the optical regime this theory was put on a 
rigorous foundation by Gardiner and CoUett [^Sl- This theory is based on 
an electric-dipole coupling, and involves key approximations which rcly on 
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a ràpid oscillation (at optical frequency ujo) of the system dipole due to the 
System Hamiltonian Hq. Let that dipole, in the interaction picture of Hq, be 
proportional to 

c{t) exp(— i^oí) + c^t) exjp{ÍLüQt). (1-87) 

Here c and are dimensionless Heisenberg-picture system operators, nor- 
malized so that cc^ |0) = |0), where |0) is the ground state of Hq. 

Let the system be located at z = and consider an incoming (z < 0) and 
outgoing (z > 0) continuum field with operator ò(z, t) as defined previously 
Then, under the rotating wave approximation 2 , the dipole coupling to the 
field at z = can be modeled by the Hamiltonian 

Vit) = z^[c{t)b\0, t) - 6(0, í)ct(í)]. (1.88) 

Here Ti = 1 and 7 ^ is the characteristic decay rate of the system. For a 
general interaction one would have to use both space and time arguments for 
the continuum field. However, the coupling considered here is strictly local. 
Both before [z < 0) and after (2 > 0) the interaction, the field still freely 
propagates. This means that, although it is necessary to use a spatial as well 
as a temporal argument, the spatial argument need only have two vàlues: 
before and after. In order to conform with pre-established usage 123', these 
vàlues will be called in and out The input field is defined (in the Heisenberg 
picture) as 

òin(í) =&(0-,í), (1.89) 
and the output field as 

bout{t) = b{0+,t). (1.90) 
Both of these fields obey the commutation relations 

m,b\t')] = 5{t-t'), (1.91) 

at least for time differences shorter than the delay time in any feedback loop. 

The interaction of the external field with the system at time t is described 
by the Hamiltonian (|1.88l) . It is necessary to be careful in dealing with this 
because of the singular commutation relations (|1.91l) . A convenient way is to 
use the quantum Itò stochastic calculus In the Heisenberg picture, the 
infinitesimal unitary transformation pertaining to the interaction at time t is 

U{t,t + dt) = exp ^c{t)dBl{t) - y/^c{t)UBin{t) - iH{t)dt , (1.92) 

where 7í is a system Hamiltonian representing any small perturbation on top 
of the free energy (such as classical driving), and 

dBin{t) ^ bin{t)dt (1.93) 



1 Squeezing and feedback 21 



is the quantum analogue of the Wiener increment Note that the field 
operator in the unitary transformation 11.92(1 is the bath immediately before 
it interacts with the system, rather than the bath at the instant it is inter- 
acting with the system. This is the appropriate operator if one treats the 
stochastic increment dBin{t) in the Itò sense. This means that (11.92(1 must 
be expanded to second order in the increment. For an input in the vacuum 
state, the operator di?in has a vanishing first order moment, and a single 
nonvanishing second-order moment Í2R| 

dBi^{t)dBl{t) = dt. (1.94) 

This could be thought of as vacuum noise. 

Applying the unitary transformation ((1.92(1 to an arbitrary system oper- 
ator s{t) yields 

s{t + dt) = U\t, t + dt)s{t)U{t, t + dt) 

= s + 7 (c^sc — jc^cs — jsc^c) dt 

- [-iHdt + ^dBl{t)c - ^cUBUt),s], (1.95) 

where all system operators in ((1.95() have the time argument t. This is what 
I wih refer to as a quantum Langevin equation (QLE) for s. Because di?in(í) 
is the bath operator before it interacts with the system, it is independent of 
the system operator s(í). Hence one can derive 

(s) = ( [7 {chc ~ ictcs - isc^c) + i[H, s]] ) . (1.96) 

From this it apparent that there is a Schròdinger picture representation 
of the system dynamics. In this picture 

(s(í)) -Tr[p(í)s], (1.97) 
where p(t) is the state matrix for the system which evidently obeys 

p = -fV[c]p-i[H,p], (1.98) 
where for arbitrary operators A and B 

V[A]B = ABA^ - \A^AB - ^BA^A. (1.99) 

An equation of the form ((1.98(1 is known as a master equation. I will discuss 
the master equation further in Sec. I1.6.Ï1 

Although the noise terms in ((1.95() do not contribute to ((1.96II . they are 
necessary in order for ((1.95(1 to be a vahd Heisenberg equation of motion. If 
they are omitted then the commutation relations of the system will not be 
preserved. As well as giving the evolution of the system, p.92(l transforms 
the input field operator into the output ficld operator: 

6out(í) = U\t, t + dt)binit)U{t, t + dt)= hn{t) + y/jc{t). (1.100) 
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Assuming that òout (t) has a large coherent amphtude then one can then deriva 
an expression for its hnearized intensity fluctuations, which are proportional 
to 

^out(í) = Xi„(í) + V^^x{t), (1.101) 

where Xo^tit) is thc fluctuation quadratura oparator for the output field as 
usual, and whara Sx — x — {x)ss, where x — c + c' is a systcni quadratura 
with maan steady-state value (a;)ss- 



1.5 Feedback based on Nonlinear Measurements 
1.5.1 QND Measurements 

Saction 11.31 showed that it was not possible to create squeezed Ught in the 
conventional sense using ordinary photodatection and hnear faedback. Al- 
though the quantum theory appeared to show that the hght which fah on 
tha detector in the faedback loop was sub-shot noise, this could not be ex- 
tractad because it was damohshed by tha detector. An obvious way around 
this would be to usa a quantum non-demolition quadrature (QND) detector. 
One way to achiave this is for two fialds of different fraquency to interact 
via a nonlinear rafractive index. In order to obtain a large nonlinearity, larga 
intensities are raquired. It is aasiest to build up large intensities by using a 
resonant cavity. To describa this requiras tha theory of quantum Langevin 
equations just presented. 

Consider tha apparatus shown in Fig. 11.21 Tha purposa of the detector 
in tha feedback loop is to maka a QND measuramant of the quadrature X^^ 
of tha fiald &in. This fialds drivas a cavity moda with dacay rata n dascribad 
by annihilation oparator a. This moda is couplad to a second mode with 
annihilation oparator c and decay rate 7. Ignoring practicalities, I will taka the 
coupling batwean tha two modes to be the ideal QND coupling, as considered 
initially by Hillery and ScuUy [31] and Yurke namely 

ií=|xV, (1-102) 
whara 

2;"=a + a^; = -ic + ic'' . (1.103) 

As described in |36| . this Hamiltonian could in principia ba realized by a 
crystal with a x^^^ nonlinearity, combined with other processes. This Hamil- 
tonian commutes with tha x"" quadratura of moda a, and causes this to driva 
tha x'^ quadratura of mode c. Thus measuring tha AT^^j quadrature of the 
output fiald dout from moda c will give a QND measurement of a + a^, which 
is approximately a QND measurement of X^^. 
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Fig. 1.2. Diagram for a traveling-wave feedback experiment based on a QND mea- 
surement. Traveling fields are denoted b and d. The first beam splitter transmittance 
rji is variable. A cavity (drawn as a ring cavity for convenience) supports two modes, 
a (sòlid line) and c (dashed line). The decay rates for these two modes are k and 
7 respectively. They are coupled by a nonlinear optical process indicated by the 
crystal labeled x- The perfect homodyne detection at the detector labeled h yields 
a photocurrent proportional to X^^^^ = dout + dl^t 



Following the analysis of |36| . the quantum Langevin equations for the 
quadrature operators are 

x"" ^-^x'^-V^X^^, (L104) 

i-^-lx^-^Xf^ + Xx". (1.105) 

These can be solved in the frequency domain as 

^ _ V«|àM (1.106) 

X (oj) = —- — . (1.107) 
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The quadrature of the output field dout = cíin + is therefore 

Xouti^) , (1.108) 

where I have defined a quality factor for the measurement 

Q = 4x/V7^. (1.109) 
In the limits w ^ k, 7 and Q 3> 1 we have 

xit(^)^-Q^iU'^), (1-110) 

which shows that a measurement (by homodyne detection) of the X quadra- 
ture of dout does indeed effect a measurement of the low frequency variation 
in XP,. 

To see that this measurement is a QND measurement we have to calculate 
the statistics of the output field from mode a, that is bo^t- From the solution 
we find 

xU-)^-^^^^xU^). (1.111) 

k/2 + loj 

That is, apart from an irrelevant phase factor the output field is identical to 
the input field, as required for a QND measurement. Of course we cannot 
expect the other quadrature to remain unaffected, because of the uncertainty 
principle. Indeed we find 

~, (^-2^^) -, Q^KY,ÍÍLü)/{^ + 2t^) 

which shows that noise has been added to Yi^. Indeed, in the good measure- 
ment limit which gave the result 1)1. 110(1 we find the phase quadrature output 
to be dominated by noise: 

tU^) ^ QY^iiu;). (1.113) 



1.5.2 QND-based Feedback 

We know wish to show how a QND measurement, such as that just considered, 
can be used to produce squeezing via feedback. The physical detalls of how 
the feedback can be achieved are as outlined in Sec. 11.3.11 In particular, in 
the limit that the transmittance ryi of the modulated beam splitter goes to 
unity, the effect of the modulation is simply to add an arbitrary signal to the 
amplitude quadrature of the controUed beam. That is, the modulated beam 
can be taken to be 



bi{t) = bo{t ^ Ti) + /3Sm{t - Ti), 



(1.114) 
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where Òq is the beam incident on the modulated beam splitter, as in Sec. ll.3.H 
In the present case, is then fed into the QND device, as shown in Fig. ll.2l 
so 

&in(í) =Òl(í), (1.115) 

and the modulation is controlled by the photocurrent from a (assumed per- 
fect) homodyne measurement of ^out- 

/>OC 

5vi{t) = ^J^ mxLtit -ro- s)ds. (1.116) 

Here tq is the delay time in the feedback loop, including the time of flight 
from the cavity for mode c to the homodyne detector, and h{s) is as before. 

Substituting (|1.114|l - 11.116|l into the results of the preceding subsection 
we find 

-^out[^ ; , : — 1 (-'-■-'--L'j 

where T = ti + tq is the total round-trip delay time, and 

represents the frequency response of the two cavities. If we assume that the 
field din is in the vacuum state then we can evaluate the spectrum of amph- 
tude fluctuations in X^^^. to be 

sLi^) = + . (1.119) 

|1 ^ 9P{^)h{ijj) exp(— ïCijr)p 

Clearly in the hmit Q ^ oo, the added noise term in the amphtude 
spectrum can be ignored. Then for sufíiciently large negative g the feedback 
will produce a sub-shot noise spectrum. Note the difference between this case 
and that of Sec. 11.31 Here the squeezed hght is not part of the feedback loop; 
it is a free beam. The ultimate limit to how squeezed the beam can be is 
caused by the noise in the measurement. In the limit 5 — )■ — cx) we find 

5„^„t(w)mi„ = \QP{u)h{u)\-\ (1.120) 

Since 6out is a free field, not part of any feedback loop, it should obey 
the Standard commutation relations. This is the case, as can be verified from 
the expression (|1.112|l for Y^^^ (which is unaffected by the feedback). Conse- 
quently, the spectrum for the phase quadrature 

SlM = Sl{^) + \QP{üj)\^ (1.121) 

shows the expected increase in noise. It is not difficult to show that the uncer- 
tainty product S^^^{íjj)S^^^{íjj) is greater than or equal to unity, as required. 
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1.5.3 Parametric Down Conversion 

The preceding section showed that feedback based on a perfect QND mea- 
surement can produce squeezing. This has never been done experimentaUy 
because of the difhcuhy of building a perfect QND measurement appara- 
tus. However, it turns out that QND measurements are not the only way to 
produce squeezing via feedback. Any mechanism which produces correlations 
between the beam of interest and another beam which are more "quantum" 
than the correlations between the two outputs of a Hnear beam sphtter can 
be the basis for producing squeezing via feedback. Such a mechanism must in- 
volve some sort of optical nonhnearity, which is why this section V is entitled 
"Feedback based on Nonhnear Measurements" . 

The production of amphtude squeezing by feeding back a quantum - cor- 
related signal was predicted [SZ| and observed [HH| by Walker and Jakeman in 
1985, using the process of parametric down-conversion. An improved feedback 
scheme for the same system was later used by Tapster, Rarity and Satchell 
| 39| to obtain an inferred amplitude spectrum {Lü)aún — 0.72 over a limited 
frequency range. 

The essential idea is as foUows. A non-degenerate optical parametric oscil- 
lator (NDOPO) can be realized by pumping a crystal with a x^^^ nonhnearity 
by a làser of frequency ui and momentum ki . For particular choices of cui the 
crystal can be aligned so that a pump photon can be transformcd into a pair 
of photons each with freqüències uj2, lu^ and momenta k2, ks with lui — ^^2+^3 
and ki = \í2 + h^. A pair of down-converted photons will thus be correlated 
in time because they are produced from a single pump photon. On a more 
macroscopic level, this means that the amplitude quadrature fluctuations in 
the two down-converted beams are positively correlated, and the coefScient of 
correlation can in principle be very close to unity. In this ideal limit, measur- 
ing the intensity of beam 2 should give a readout identical to that obtained 
from measuring beam 3. In effect, the measurement of beam 2 is like a QND 
measurement of beam 3. Thus, feeding back the photocurrent from beam 2 
with an negative gain should be able to reduce the noise in beam 3 below the 
shot-noise limit. In the limit of perfect detection, the optimum gain becomes 
arbitrarily large. 

In |37l39j the negative feedback was effected by controUing the power of 
the pump làser (which controls the rate at which photon pairs are produced) . 
This maintains the symmetry of the experiment so that the fed-back pho- 
tocurrent from beam 2 has the same statistics as the photocurrent from the 
free beam 3. In other words, they will both be below the shot-noise, whereas 
without feedback they are both at best shot-noise limited. However, it is not 
necessary to preserve the symmetry in this way. The measured photocur- 
rent from beam 2 can be fed forward to control the amplitude fluctuations 
in beam 3 (for example by using an electro-optic modulator as described in 
Sec. I1.3.ï|l . This feedforward was realized experimentaUy by Mertz et al. in 
1990 ||4(}j achieving similar results to that obtained by feedback The les- 
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son is that unless one is concerned with light inside a feedback loop, there is 
no essential difFerence between feedback and feedforward. Indeed, the squeez- 
ing produced by QND-based feedback discussed in Sec. 11.5.21 could equally 
well have been produced by QND-based feedforward. 

1.5.4 Second Harmònic Generation 

A third example of using quantum correlations between two beams to good 
efïect in a feedback (or feedforward) loop is that of second-harmonic gen- 
eration. Like parametric down conversion, this involves a x^^-* crystal, but 
here the two low frequency modes are assumed to be degenerate, so that 

= ix>z and k2 = ka, and the polarizations are the same also. For simplic- 
ity, I will call the low frequency mode the red mode and the high frequency 
mode with frequency íü\ — 2uj2 the green mode. For second harmònic gener- 
ation, it is the red mode which is pumped (that is, the reverse of parametric 
down-conversion) . In this case, the pumped mode (red) is also treated as a 
quantum system so that the two beams are the red and green beams. The 
system is easily treated using quantum Langevin equations 

The proposal made in 0^ is as foUows. The first harmònic (red) resides 
in a good cavity and is driven by a làser at the end with very high reflectivity. 
The second harmònic (green), which is generated in the red cavity, is refiected 
at one end only (giving two passes) so that it forms a single output beam. 
By itself (that is, in the absence of feedback) this system produces amplitude 
squeezing in both the first and second harmònic |44| . The amplitude squeez- 
ing in the red can be understood to be due to two-photon absorption, which 
is the effect of the adiabatically eliminated green mode. Such nonlinear ab- 
sorption preferentially damps large fluctuations and so reduces the variance. 
The antibunching in the green mode can be attributed to the fact that the 
creation of a green photon requires the loss of two red photons, which reduces 
the chances of this event re-occurring. 

Without feedback, the optimum low frequency squeezing in the red mode 
output is 



Here x is the ratio of the nonlinear intensity loss rate to the linear intensity 
loss rate for the red mode. Now there are two possible ways to feed back 
onto the driving of the system. The first is to use the green photocurrent to 
control the amplitude of the red driving, in order to try to reduce the noise 
in the red output light. A linearized treatment gives a new minimum of 



(1.122) 



and that in the green mode is 



Sg{0)=g for X^^- 



(1.123) 



^.''(0) - - for X = 1, 9 



1 

2' 



(1.124) 
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which is significantly less than the no-feedback value 122(1 . Here, g is the 
low frequency loop transfer function, equal to the round-loop gain of the 
feedback Note that it is positive, corresponding to destabihzing feedback. 
This is difïerent from the cases of QND-based feedback and down-conversion- 
based feedback. 

The other sort of feedback is to control the driving of the red mode using 
the detected red photocurrent, looking to enhancc the squeezing in the green 
output. It turns out that the minimum 5*^(0) of H1.123|l as x ^ co cannot 
be lowered by feedback. However, for finite vàlues of x, the feedback can give 
a definite improvement. For example, with no feedback 

Sf{0) = i for X - 1, (1.125) 
whereas the optimal feedback gives 

5f (0) = i for X = 1, 9^\- (1.126) 

All of these results are calculated for the case of unit efficiency photode- 
tectors. The effect of non-unit efhciency is to reduce the effectiveness of the 
feedback, and to alter the conditions of optimality. The general solution, 
including the full spectrum with an arbitrary transfer function gh{üj), is con- 
tained in 



1.6 Quantum Trajectòries 
1.6.1 The Master Equation 

The Quantum Langcvin Equation (QLE) considered in Sec. 11.41 are Heisen- 
berg picture equations which detall the effect of the bath on the system. They 
also, through the input-output relations, specify the effect of the system on 
the bath, and hence the relation between the system evolution and the mea- 
sured photocurrents. For systems which have linear QLEs, this is the easiest 
method to analyze their evolution. However, even some simple systems (such 
as a two-level atom) do not have linear QLEs. Also, sometimes a better in- 
tuition about feedback can be gained by working in the Schròdinger picture 
rather than the Heisenberg picture. 

If one is interested only in the evolution of the system then there is a sim- 
ple Schròdinger picture equivalent to the QLE. This is the quantum master 
equation, an example of which was already derived as 1)1. 98|) from the corre- 
sponding QLE. Unlike a QLE, a master equation is always a linear equation. 
It is possible to derive the master equation directly from the system-bath cou- 
pling without deriving the QLE, simply by tracing over the bath. I present 
this derivation below, because it shows clearly that this step (ignoring the 
bath) is not an essential one. If instead of ignoring the bath, one measures 
it, then onc obtains a quantum trajectory equation. This is a different sort 
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of Schròdinger picture equivalent to the QLE which is more general than 
the master equation, as it can also relate the system evolution to the mea- 
sured photocurrents. This is necessary if one is to consider feeding back the 
photocurrent to alter the dynamics of the system. 

For convenience, I will work in the interaction picture, rather than the 
Schròdinger picture, so that the free evolution causing the oscillation of the 
dipole can be ignored. The interaction between the system and the bath is 
now given by 

Vit) = ihy/^ [b^t)c - c^bit)], (1.127) 

where c is a slowly-varying interaction picture operator, and b{t) is also a 
slowly-varying interaction picture bath operator. The timc-depcndence is 
maintained for the bath operator, because the free Hamiltonian of the bath 
causes propagation at the speed of light, so a new part of the bath interacts 
with the system at each new point in time. These parts are labeled by the 
time of interaction t. Thus, b(t) is an operator in the Hilbert space for a 
particular part of the bath. Each part has its own state matrix For the 
incoming field to be a bath requires that its total state matrix be the direct 
product of the state matrices of the parts . That is to say, the temporally 
separate parts of the bath must be unentangled. 

Let the system at time t be known to be p(t). Thus the initial state of 
the system and (relevant part of the) bath at time t is 

R{t) ^ ^i{t)® p(t). (1.128) 
The infinitesimally evolved state is 

R{t + dt) = U{t,t + dt)[fi{t) ^ p{t)]U\t,t + dt), (1.129) 
where 

U{t, t + dt) — exp [y^dB\t)c- c^y^dB{t)-iHdt\ , (1.130) 

where H is the system Hamiltonian in the interaction picture. If the input 
bath is in the vacuum state then all of the first and second order moments 
of dB{t) vanish except, as in H1.94|l . 

Tr[dB{t)dB''{t)p{t)] = dt. (1.131) 

Thus, it is necessary to expand some of the terms in U{t,t + dt) to second 
order. The result for R{t + dt) is 

m(í) «> P(í) + V7 [dB\t)c-cUB{t),fi{t)® p{t)] (1.132) 
+-i {dB\t)n{t)dB{t) ® cp{t)c^ - ^dB{t)dB\t)p{t) ® c^cp 
- \p{t)dB{t)dB\t) ® pc^ c] -i[Hdt,p{t) (g> p{t)]. 
The infinitesimally evolved reduced state matrix for the system is given by 
p{t + dt) = Tr^[i?(í + dt)]. (1.133) 
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Taking the trace over the bath state in p.l32|l yields the master equation 

p = ^V[c]p-i[H,p]. (1.134) 

The first term here represents irreversible evolution (damping in this case), 
and is of the unique form for such evolution as proved by Lindblad '45' . The 
second term represents unitary evolution which is of course reversible. 



1.6.2 Photon Counting 

To consider measuring the output bath it is useful to explicitly write down its 
state. Before interacting with the system, the input bath state is p{t) — |0) (0| 
for all time í, where |0) is the lowest eigenstate for a^a. Here, a = y/dtbo{t) — 
dBo{t)/Vdt, so that a^a has integer eigenvalues representing the number of 
photons arriving in the interval of time [t, t+dt). Let the state of the system at 
time t be p(t), independent of The entangled state after the interaction 
of duration dt is, from (|1. 132(1 

R{t + dt) ^ |0)(0| (g)p(í) + V7^[|l>(0| «)cp(í) + \Q){l\(g)p{t)c^ 
-idt\0){0\ ® [H,p{t)] + -fdt{\l)(l\(i)cp{t)c'' 
~ i|0)(0|(» [c^cp{t)+pit)c^c]}. (1.135) 

The free dynamics of the electromagnètic field will now remove the bath state 
from the system, so that the entanglement in H1.135|l will be maintained. If 
the outgoing bath is ignored, then the system obeys the master equation 

To obtain Information about the system, the outgoing field must be mea- 
sured. The obvious measurement to consider is photon counting. This I will 
model by projecting the field states onto eigenstates of a^a. It is possible to 
consider specific models for photon detectors, usually based on atòmic Sys- 
tems [21 . However, these simply remove the measurement step, where possibil- 
ities become actualities, one step further along the von Neumann chain |46| . 
There is little which results from such models which cannot be achieved by in- 
cluding losses and convolving the classical photocurrent with an empirically 
derived detector response function. Therefore, I will simply use projection 
measurement operators 

Po = |0)(0|; Pi = |l>(l|. (1.136) 

If there is a nuU count then the unnormalized state matrix of the system 
and bath (whose norm gives the probability of this outcome) is given by 

Ro{t + dt) = PoR{t + dt)PQ ^ Po(g) poit + dt), (1.137) 

where 



po{t + dt) = p{t) - 7Í{c^c, p{t)}dt - i[H, p{t)]dt. 



(1.138) 
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This has a norm only infinitesimally difïerent from one, so for almost all time 
intervals no photons are detected in the output. 

If a photon is detected, then Ri{t + dt) = Pi (E) pi{t + dt). That is, the 
systeni jumps into the unnormalized conditioned state 

pi{t + dt) = -fcp{t)cUt. (1.139) 

The norm of this state matrix is equal to the probability of a detection 
occurring in the interval [t,t + dt), and is equal to ^{c^c)dt. Clearly, the 
unconditioned master equation evolution H1.134|l is retrieved by averaging 
ovcr the two possible results 

p{t + dt) = po{t + dt) + + dt). (1.140) 

It is useful and elegant to reformulate this evolution in the form of an 
explicitly stochastic evolution equation. This equation of motion specifies the 
quantum trajectory of the system. Since the measurement result is a point 
process, it can be represented by a random variable dNc{t) representing the 
increment (either zero or one) in the photon count in the interval [t,t + dt). 
It is formally defined by 

E[dN,{t)] = Tr[c^cp,{t)]-fdt, (1.141) 
dNc{t)^ = dNcit). (1.142) 

Here the subscript c indicates that the quantity to which it is attached is 
conditioned on previous measurement results, arbitrarily far back in time. 
The conditioned state matrix obeys the stochastic master equation (SME) 

dp^{t) ^ {dNc{t)g[c] ~ dtn [íH + Í7ctc] } pc(í)· (1.143) 

Here, the nonlinear (in p) superoperators Ç and TL are defined by 

^Hp^-^P^-P, (1.144) 

H[r\p = rp + pr^ - Tr[rp + pr'^]p. (1.145) 

The nonlinearity of the SME H1.143I) is indicative of the fundamental non- 
linearity of quantum mcasurements. The original master equation for p[t) = 
E[pc{t)] can be restorcd simply by replacing dNc{t) in (|1.143|1 by its cnscmble 
average value H1.141|l . 

Because of the assumed perfect detection, the stochastic equation for the 
state matrix is equivalent to a stochastic equation for the state vector. The 
unraveling of the master equation into a stochastic Schròdinger equation 
representing photon counting is the most commonly used quantum trajec- 
tory for numerical simulations 47 48 49 50 51 52 . From the point of view of 
measurement theory, the stochastic master equation is of more use than the 
stochastic Schròdinger equation because it is more transparently related to 
the unconditioned master equation and because it can be generalized to cope 
with ineSicient detectors. This will be dealt with explicitly for the case of 
homodyne detection. 
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1.6.3 Homodyne Detection Theory 

The unraveling of the master equation (11.13411 as a quantum trajectory is not 
unique. DifFerent detection schemes will result in difFerent quantum trajectory 
equations. For squeezing, the most useful detection technique is homodyne 
detection. In the simplest configuration, the output field of the cavity, &out = 
V + y/^c, is sent through a beam spHtter of transmittance 77 very close to 
one. Into the other input port of the beam sphtter is injected a very strong 
coherent field. This has the same frequency as the system dipole, and is 
known as the local oscillator. The transmitted field is then represented by 
the operator 

bi = v + ^c + P, (1.146) 

where /? is a complex number representing a coherent amplitude, such that 
|/3p/(l — 77) is equal to the input photon flux of the local oscillator. The 
photodetection operators are then applied as above, with the annihilation 
operator defined as a = hi^fdt. 

Let the coherent field /3 be real, so that the homodyne detection leads to 
a measurement of the x quadrature of the system dipole. Also, let us measure 
time in units of 7^^ so that this parameter disappears from our equations. 
Then the rate of photodetections at the (perfect) detector 

E[d7Vc(í)] = Tr[(/32 +I3x + c'^ c)pc{t)]dt. (1.147) 

where x — c-\- as previously. In the limit that /3 is much larger than c, this 
rate consists of a large constant term plus a term proportional to plus a 
small term. It is not difficult to show that the SME for the conditioned state 
matrix is altered from (|1.143|l to 

dp^{t) = {dNS)g[c + 13]+ dtn [-~iH - (3c- \c^c\ } pc(í)- (1.148) 

The ideal limit of homodyne detection is when the local oscillator ampli- 
tude goes to infinity. In this limit, the rate of photodetections goes to infinity, 
but the effect of each on the system goes to zero, because the field being de- 
tected is almost entirely due to the local oscillator. Thus, it should be possible 
to approximate the photocurrent by a continuous function of time, and also 
to derive a smooth evolution equation for the system. This was done first by 
Carmichael ^Sl- A more rigorous working of the derivation is found in j5;-i| . 
The result is that over a time much longer than but much smaller 

than unity, the system evolution can be approximated by the SME 



dpc{t) = -i[H, p^{t)]dt + V[c\pc{t)dt + dW{t)n[c\pc{t). (1.149) 

Here dW{t) is an infinitesimal Wiener increment 54 satisfying 

E[dW(t)] = 0, (1.150) 

dW{t)^ = dt. (1.151) 
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Thus, the jump evolution of l|1.148|l has been replaced by difFusive evolution. 
Equation 1)1.149(1 is, by its derivation, an Itò stochastic master equation, 
where the equal-time stochastic increment dW{t) is independent of the state 
of the system Pc{t). It is trivial to see that the ensemble average evolution 
reproduces the nonselective master equation (|1.134l) by eliminating the noise 
term. 

Just as the /? — > oo leads to continuous evolution for the state, it also 
changes the point process photocount into a continuous photocurrent with 
white noise. Removing the constant local-oscillator contribution and scaling 
appropriately gives 

/ho-(í) ^ hm lim ^^^W -/^^'^^ = + ^(,), (1.152) 

ít— >0/3^oo pot 

where Ç(í) = dW{t)/dt. It is not difficult to see that if the detector efficiency 
is 77, the homodyne photocurrent becomes 

/^-(í) ^ hm hm ^ Vv{^)c{t)+m, (1.153) 

(5í^o/3^oo y/rj pot 

and the SME (|1.149|l is modified to 

dp = -i[H, pc{t)]dt + V[c]pc{t)dt + y/7jdW{t)n[c]pc{t). (1.154) 



1.6.4 Homodyne-mediated Feedback 

In the foUowing section I will consider the use of feedback to produce or 
enhance squeezing in intracavity fields. Since squeezing is the reduction in 
fluctuations of one quadrature of a field, the obvious sort of feedback to 
consider is one using the homodyne photocurrent obtained from measuring 
the output field. Here I will develop the theory for describing this sort of 
feedback. As well as being more relevant for our purposes than feedback 
using the direct detection photocurrent jEli it is also somewhat easier to 
treat theoretically, which is why it was derived first |l()lll| . 

In principle, the homodyne photocurrent could be subject to any sort 
of filtering prior to being fed back, including nonlinear filtering. However it 
turns out that for the applications we wish to consider, only linear filtering 
is desired. Also, rather than using a response function h{t) I will simply take 
the feedback to be delayed by a time T. That means that the evolution due 
to the feedback can simply be written as 

[Pc(í)]fb = /^'"(í - T)ICpAt)/Vv, (1-155) 

where /C is a superoperator. Since /^°™(í) may be negative, /C must be such 
as to give vàlid evolution irrespective of the sign of time. That is to say, it 
must give reversible evolution with 



JCp=-i[F,p] 



(1.156) 
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for some Hermitian operator F. In other words, we can represent the efFect 
of the feedback by the Hamiltonian 

Híi,^F [{c + c))S -T) + e(í - T)/^] . (1.157) 

Because the stochasticity in the measurement (|1. 15411 and the feedback 
p.l55() is Gaussian white noise, it is relatively simple to determine the efFect 
of the feedback. Bearing in mind that the feedback must act after the mea- 
surement, and that 1)1.155(1 must be interpreted as a Stratonovich equation 
|11| . the resuh for the total conditioned evolution of the system is 

pM + dt) = \ 1 + ICUc + c^)M - T)dt + dW(t - T)/J7Ï] + —IC^dt 

X {1 + n[-iH]dt + V[c]dt + y^dW{t)n[c\}pcit). (1.158) 
For T finite, this becomes 

1 



dp^it) = dt ^n[-iH] + V[c] + {c + cT)c(< - T))C + 2^^^| Pc(í) 

+ dW{t - T)ICpc{t)/^ + ^dW{t)n[c\pc{t). (1.159) 

On the other hand, putting T = in (|1.158|l gives 

dpc{t) = dt {-i[H, pc(í)] + 25[c]/5c(í) - i[F, cpc{t) + Pc{t)c''] 

+ V[F]p4t)/Tj} + dW{t)H[^c ~ iF/^]p,{t). (1.160) 



1.7 Intracavity Squeezing 
1.7.1 The Linear System 

In order to understand the effect of quantum limited feedback on intracavity 
squeezing, it is useful to consider an exactly solvable system. In this section, I 
will mainly be foUowing jll| in considering the case of a linear optical system, 
with linear feedback based on homodyne (or QND) detection. By a linear sys- 
tem, I mean that the equation of motion for the two quadrature operators 
are linear. This is approximately the case for many quantum optical systems, 
in the limit of large photon numbers. For specificity, I will chose a system 
which is exactly linear. If, as in the remainder of this chapter, one is inter- 
ested in the behaviour of one quadrature only (here the x quadrature), then 
all linear dynamics can be composed of damping, driving, and parametric 
driving. Damping will be assumed to be always present (as necessary to do 
feedback or obtain an output from the cavity) and will have rate 1. Constant 
linear driving simply shifts the origin away from a; = 0, and will be ignored. 
Stochastic linear driving in the white noise approximation causes diffusion 
in the X quadrature, at a rate l. Finally, if the strength of the parametric 
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driving {H xy) is 9 (where 0=1 would represent a degenerate parametric 
oscillator at threshold), then the master equation for the system is 

p = V[a]p + \lV[a^ - a]p + \0[a'^ - a^^, p] = CqP, (1-161) 

where a is the annihilation operator for the cavity mode. 

An alternative definition for the Unearity of the x quadrature dynamics 

is that the marginal distribution of the Wigncr function for x (which is the 
true probabihty distribution for x) obeys an Ornstein-Uhlenbeck equation. 
That is to say, 

P(x) = {d^kx + \Ddl) P{x), (1.162) 

where k and D are constants. The solution of this equation is a Gaussian 
with variance 

y = ^. (1.163) 

For the particular master equation above (the properties of which will be 
denoted by the subscript 0), the drift and diffusion constants are 

fco = i(l + 0), (1.164) 
Do = 1+1. (1.165) 

In this case, Vq = (1 + Z)/(l + 6*). If this is Icss than unity, the system exhibits 
squeezing of x. It is more useful to work with the normally ordered variance, 
which becomes negative if the x quadrature is squeezed. Here, I will denote 
it 

U = V-l, (1.166) 
which for this system takes the value 

í/o=[^. (1.167) 

If the system is to stay below threshold (so that the y quadrature does not 
become unboundcd), thcn the maximum value for 6 is onc. At this value, 
Uç) = —1/2 when the x diffusion rate l = 0. Therefore the minimum value of 
squeezing which this linear system can attain as a stationary value is half of 
the theoretical minimum oiU = — 1. 



1.7.2 Homodyne-Mediated Feedback 



We now wish to consider the cffcct of homodyne-mediatcd feedback on the 
intracavity light. This is most easily understood using the quantum trajectory 
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picture in the Markovian limit. Thus we want the stochastic master equation 
for the conditioned state matrix Pc{t) 11.160(1 

dpc{t) = dt (^CoPc{t) +IC[apc{t) + Pc{t)a''] + ^/C^pc(í) 

+ dWit){^n[a]+)C/^)pcit). (1.168) 

Here, £o is as defined in 1(1.161(1 . 

The question now arises as to what to choose for the JC. Seeking to reduce 
the fluctuations in x suggests the feedback operator, which is related to /C by 
1(1.156(1 . should be 

F = -Xy/2. (1.169) 

As a separate Hamiltonian, this translates a state in the negative x direction 
for A positive. By controlHng this Hamiltonian by the homodyne photocurrent 
one thus has the abihty to change the statistics for x and perhaps achieve 
better squeezing. This Hamiltonian can be effected by driving the cavity (at 
a second mirror which can be assumed to have a negligible loss rate compared 
to the first mirror). Using this choice and changing ((1.168(1 into a stochastic 
Liouville equation for the conditioned Wigner function gives 



dPc{x) ^ dt 



d^iko + X)x + (Do + 2X + Xyrj) 



Pc{x) 

+ dW{t) [07 (x - x,{t) + d^) + (A/V^)a,] P,{x), (1.170) 

where íc(í) is the mean of the distribution Pc{x) and dW{t) is as usual. 

This equation is obviously no longer a simple Ornstein-Uhlenbeck equa- 
tion. Nevertheless, it still has a Gaussian as an exact solution, as can be 
shown by direct substitution. The mean Xc and variance of the condi- 
tioned Gaussian distribution are found to obey 

ée = -(fco + \)x,+m [^{Vc - 1) - (A/V^)] , (1.171) 
V, = -2kç,V, + D^-T^{V,~\f . (1.172) 

Two points about the evolution equation for Vc are worth noting. It is com- 
pletely deterministic (no noise terms) , and it is not influenced by the presence 
of feedback. Furthcrmore, for this linear system, it is independent of Xc- Thus, 
the stochasticity and feedback terms in the equation for the mean do not even 
enter that for the variance indirectly. 

The equation for the conditioned variance is more simply written in terms 
of the conditioned normally ordered variance Uc = Vc — 1 

Úc = -2fcoí7c - 2ko + Do~ rp^. (1.173) 

On a time scale as short as a cavity lifetime, Uc will approach its stable 
steady-state value of 



(1.174) 
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Substituting the steady-state conditioned variance into 1)1.17111 gives 



ïí one were to choose 



-ka + ^ kl + i-i{-2ka + Da) - X 



X 



-ka + \l kl + T^{-2ka + Da) 



-ka 



ka 



2r]koUo 



(1.175) 



(1.176) 



then there would be no noise at all in the conditioned mean and so one could 
set Xc = 0. In other words, this value of A is precisely the value required 
to minimize the unconditioned variance under feedback. When all fluctua- 
tions in the mean are suppressed, the unconditioned variance is equal to the 
conditioned variance. 

In general, the unconditioned variance will consist of two terms, the con- 
ditioned quantum variance in x plus the classical (ensemble) average variance 
in the conditioned mean of x: 



Ux^U, + E[xl]. 
The latter term is found from H1.175|l to be 



E[xl 



' 2(fco + A) 
Adding (|1.174|l gives 



-(fco + A) + ^ kl + r]{-2ka + Da) 



Ux = V 



_^X^ +f]{-2ko + Pa) 
2{ko + X) 



(ka + Xy 
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(1.177) 



(1.178) 



(1.179) 



An immediate consequence of this expression is that U\ can only be neg- 
ative if Uq is. That is to say, classical feedback based on homodyne detection 
cannot produce intracavity squeezing. However, this does not mean that the 
feedback cannot enhance squeezing. Obviously, the best intracavity squeezing 
will be when 77 = 1, in which casc the intracavity squeezing can be simply 
expressed as 



Un 



ka -l 



v/1 + Ro) 



(1.180) 



where Rq = (-2fco + £'o)/fco > -1- I* can be proven that that t/min < Uo-, 
with equality only if 77 = or C/q = 0. This result implies that the intracavity 
variance in x can always be reduced by classical homodyne-mediated feed- 
back, unless it is at the classical minimum. In particular, intracavity squeez- 
ing can always be enhanced. For the parametric oscillator defined originally 
in (|1. 16111 . with l — 0, í/min = S/rj. For r; = 1, the (symmetrically ordered) 
X variance is Vmin = 1 — 0. The y variance, which is unaffected by feedback, 
is seen from H1.161|l to be (1 — é*)^^. Thus, with perfect detection, it is pos- 
sible to produce a minimum unccrtainty squeezed state with arbitrarily high 
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squeezing as 1. This is not unexpected as a parametric amplifier (in 

an undamped cavity) also produces minimum uncertainty squeezed states. 
The feedback removes the noise which was added by the damping which is 
necessary to do the measurement used in the feedback. 

The reason that this feedback cannot produce squeezing is that the con- 
ditioning of the variance according to H1.173|l cannot change the sign of the 
normaUy-ordered variance U. The homodyne measurement does reduces the 
conditioned variance, except when it is equal to the classical minimum of 1. 
The more efhcient the measurement, the greater the reduction. Ordinarily, 
this reduced variance is not evident because the measurement gives a random 
shift to the conditional mean of x, with the randomness arising from the shot 
noise of the photocurrent. By appropriately fecding back this photocurrent, it 
is possible to precisely counteract this shift and thus observe the conditioned 
variance. 

If the time delay T in the feedback loop is not neghgible then the counter- 
action wih be less than perfect. It is possible to calculate this efFect exactly 
for an arbitrary linear feedback response using the quantum trajectory theory 
However, it would generally be easier to return to the approach based on 
quantum Langevin equations 55 . For short delay T <^ 1, there is a simple 
expression for the modified normally ordered variance: 

Ux-T^Uxil + XT). (1.181) 

For squeezed systems, with Uq < 0, the optimum value of U\ occurs for A 
negative, as shown above. Thus, the time delay reduces the total squeezing 
by the factor (1 + AT). On the other hand, classical noise is reduced to í7a > 
with A positive, so that the total noise is increased by the factor (1 + AT). 
Overall, the time delay degrades the effectiveness of the feedback, as expected. 

Note that the optimal A of H1.176|l has the same sign as Uq. That is to 
say, if the system produces squeezed light, then the best way to enhance 
the squeezing is to add a force which displaces the state in the direction 
of the difference between the measured photocurrent and the desired mean 
photocurrent. This is the opposite of what would be expected classically, 
and can be attributed to the effect of homodyne measurement on squeezed 
states. For classical statistics {U > 0) , a higher than average photocurrent 
reading [^(í) > 0] leads to the conditioned mean Xc increasing (except if 
C/ = in which case the measurement has no effect). However, for nonclassical 
states with Í7 < 0, the classical intuition fails as a positive photocurrent 
fluctuation causes Xc to decrease. This explains the counterintuitive negative 
value of A required in squeezed systems, which naively would be thought to 
destabilize the system and increase fluctuations. The value of the positive 
feedback required (|1.17t)|l is such that the overall restoring force fco + is still 
positive. 

Succinctly, one can state that conditioning can be made practical by feed- 
back. The intracavity noise reduction produced by classical feedback can be 
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precisely as good as that produced by conditioning. This reinforces the sim- 
ple explanation as to why homodyne-mediated classical feedback cannot pro- 
duce nonclassical states: because homodyne detection cannot. Nonclassical 
feedback (such as using the photocurrent to influence nonlinear intracavity 
elements) may produce nonclassical states, but such elements can produce 
nonclassical states without feedback, so this is hardly surprising. In order to 
produce nonclassical states by classical feedback, it would be necessary to 
have a nonclassical measurement scheme. That is to say, one which docs not 
rely on measurement of the extracavity light to procure Information about 
the intracavity state. Intracavity measurements (in particular, quantum non- 
demolition measurements) are not limited by the random process of damping 
to the external continuum. The extra term which the measurement intro- 
duces into the nonselective master equation will not produce nonclassical 
states, but may allow the measurement to produce nonclassical conditioned 
states. One would thus expect that intracavity QND measurements would 
enable feedback to overcomes the classical limit, and I will now show that 
this is indeed the case. 



1.7.3 QND-Mediated Feedback 

The natural choice of quantum non-demolition variable is the quadrature to 
be squeezed, say x as before. I use the same model for a QND measurement as 
in Sec. II. 5. 11 Mode a is coupled to mode c by the Hamiltonian in l|l. 102(1 . The 
other dynamics of mode a are defined as before by its Liouville superoperator 
Cq. The density operator for both modes thus obeys the foUowing master 
equation: 

R = CqR - ^[x{c - c'<),R] + jV[c]R. (1.182) 

In order to treat mode c as part of the apparatus rather than part of 
the system, it is necessary to eliminate its dynamics. This can be done by 
assuming that it is heavily damped, with 7 much larger than all other rates. 
Then, apart from initial transients, it will have few photons and will be slaved 
to mode a. FoUowing Standard techniques for adiabatic elimination [^Sl gives 
the master equation for p, the density operator for mode a alone, as 

CQp + rV[x/2]p, (1.183) 

where the measurement strength parameter is _r = 4x^/7. 

Now add homodyne measurement of the b mode with efíiciency ?/. Starting 
from the conditioned state matrix Rc before the adiabatic elimination and 
foUowing it through gives the conditioning master equation for pc |11| 

dpc = dtCoPc + dtrV[x/2]pc + yGïrdW{t)n[x/2]pc. (1.184) 

Normalizing the homodyne photocurrent so that the noise is the same as in 
preceding sections gives 

Ic{t) = Vïï{x)c{t)+^{t). (1.185) 
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Here I am using H (a capital rj) for rjF as the efïective efïiciency of the 
measurement. This is related to the parameter Q in Sec. ll.5.ïl bv H = rjQ^ /A. 
Note that this is not bounded above by unity, since it is possible for to 
be much greater than onc even with x much less than 7. Recali that all rates 
are measured in units of the a mode linewidth (which was k in Sec. I1.5.H) . 

The photocurrent 111.1851) can be used in feedback onto the a mode just 
as in preceding sections. A feedback term of the form 



[Pc]íh = /c(í - T)JCpc/V}ï 

gives, in the limit T 0, the conditioned evolution 



(1.186) 



dpc = dt yCoPc + rV[x/2\pc + JCj[xpc + Pcx] 

+ dw{t) (^/mí[x/2\ + a:/\/h) pc. 



2H 



(1.187) 



Using the same expressions as in Sec. 11.7.21 implies that the probability 
distribution for the x quadrature obeys 



dP,{x) 



dt 



5,(fco + A)x+ia2(i5o + AVH) 



^c(x) 



+ dW{t) ^/YÍ[x-x^{t)] + {X/y/H)d^ P^{x) 



The mean and variance of this conditioned distribution obey 

xc = -{h + A)ïc + e(í) (Vl·lv^c - a/Vh) , 



Vc 



-2kaVc + Do- HK^ 



(1.188) 

(1.189) 
(1.190) 



These equations are identical to the corresponding equations for homodyne 
mediated feedback (|1. 171(1 and H1.172|l apart from the replacement of (14 — 1) 
by Vc and 77 by H in the measurement terms. In the limit H — > 00, H1.190|l 
predicts an arbitrarily small steady-state conditioned variance. This is char- 
acteristic of a good QND measurement. Of course, the quantum noise has 
not been eliminated but rather redistributed. For H to be large requires F to 
be large also, so that the variance in the unsqueezed quadrature is greatly in- 
creased by the measurement term in (I1.183|l . This ensures that Heisenberg's 
uncertainty principle is not violated. 

For this QND measurement the stationary value for V from 1)1.190(1 is 



(1.191) 



Thus choosing the feedback strength to be 



A = -fcn 



(1.192) 
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eliminates the stochastic element in (I1.189|l . In this case, the stationary con- 
ditioned variance p.l91|l is the minimum achievable variance. In the hmit 
H ^ oo, it is easy to see that Vmin — Ve approaches the theoretical minimum 
value of 0. That is, perfect squeezing can be produced inside the cavity by 
QND mediated feedback. In this Hmit, one requires the feedback to be very 
strong, with A ~ \/}íDq. Unhke the homodyne mediated feedback case, A 
should always be positive, as in accord with classical intuition. Indeed, all 
of the features of QND mediated feedback conform to a classical theory of 
feedback with measurements of finite accuracy (related to H) . The quantum 
nature of the feedback is manifest only in the increased fluctuations in y due 
to the measurement back-action not present classically. 

1.7.4 Mimicking a Squeezed Bath 

The application of feedback based on a QND homodyne measurement to a 
cavity State has also been considered by Tombesi and Vitali 56 . However, 
rather than directly trying to minimize the x quadrature variance of the 
field mode, as I have discussed above, their goal was to mimic the dynamics 
produced by shining a broad-band squeezed vacuum onto the cavity mirror. 
Being broad-band compared to the cavity mode, a squeezed vacuum input is 
parametrized by two numbers, M which change the single nonzero second- 
order moments of (|1.94() to 

dBl^dBin = Ndt ; dBindBl^ = {N + l)dt ; 

dBi^dBin = (dBldBl)* = Mdt. (1.193) 

Positivity of the bath state matrix requires |Mp < N{N + 1) The master 
equation resulting from such bath correlations is |2| 

p = -i[H, p] + {N + l)V[a] + NV[a^] 

-A/i[at, [a\p]] ~ Aí*i[a, [a,p]]. (1.194) 

Tombesi and Vitali show that an equation of this form can be produced using 
feedback based on a QND homodyne measurement and that, not surprisingly, 
it can can produce intracavity squeezing. 

1.7.5 The Micromaser 

A final example of feedback onto an intracavity state which can give non- 
classical noise reduction is that of the micromaser jïï7|. This consists of a 
small microwave cavity through which a monochromatic beam of resonant 
two-level àtoms is passed. The atòmic state upon exit can be measured and 
the result used in feedback. The case of modifying the cavity quality factor 
has been considered by Liebman and Milburn j58| . Because of the nature 
of the Jaynes-Cummings Hamiltonian, the micromaser dynamics are com- 
plicated without feedback, and even more complicated with. However, one 
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result is easy to explain. In the limit of short trànsit time the àtoms (as- 
sumed all to enter in the upper state) act simply as a linear amplifier of the 
cavity mode. In the absence of feedback the stationary state is thermal, with 
a photon numbcr variance much greater than the mean. With weak feedback, 
increasing the cavity damping rate whenever an outgoing atom is detected 
in the lower state, the photon distribution can be made sub-Poissonian, with 
a variance equal to half the mean. For longer trànsit times, the no- feedback 
dynamics show the effects of trapping states (where the atom undergoes an 
integer number of Rabi cycles in trànsit |S7]), and the minimum stationary 
variance is typically as low as one quarter of the mean. Feedback can pro- 
duce an arbitrary small minimum variance near a trapping state. However, 
this result is very sensitive to the trànsit time and so may be washed out by 
a realistic atòmic velocity profile. 

1.8 Feedback Master Equation 

In (ILlSUfl wc have the stochastic master equation for a system undergoing 
homodyne measurement, with instantaneous linear feedback of the homodyne 
photocurrent. This is a Markovian Itò stochastic equation, so that is possible 
to take the ensemble average simply by removing the stochastic term. This 
removes all nonlinear terms from the stochastic master equation, and gives 
the deterministic master equation 

p = -i[H, p] + V[c]p - i[F, cp + pc^] + -V[F]p. (1.195) 

V 

I will call this the homodyne feedback master equation. The first feedback 
term, linear in F, is the desired effect of the feedback which would dominate 
in the classical regime. The second feedback term causes diffusion in the 
variable conjugate to F. It can be attributed to the inevitable introduction 
of noise by the measurement step in the quantum-limited feedback loop. The 
lower the efficiency, the more noise introduced. 

The homodyne feedback master equation can be rewritten in the Lindblad 
form 0ïï| as 

p^-i\H+ i(c^F + Fc).p] + V\c - iF]p + i^PÍFlp = Cp. (1.196) 

In this arrangement, the effect of the feedback is seen to replace c by c — iF, 
and to add an extra term to the Hamiltonian, plus an extra diffusion term 
which vanishes for perfect detection. It is possible to derive an analogous 
feedback master equation for direct detection, but this is not needed for 
squeezing. In the limit where intensity fluctuations can be linearized as am- 
plitude quadrature fluctuations, direct detection is essentially equivalent to 
homodyne detection of the amplitude quadrature. 

It is very important to note that although (|1.19t)|l has the appearance of 
a normal master equation, one cannot simply use it in the customary way for 
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calculating the spectrum of the photocurrent used in the feedback loop. For 
example the two-time correlation function of the in-loop photocurrent is not 
given by the Standard expression: 

E[/h°-(í + r)l!^°"\t)] ^ íy(: x(t + T)x(t) :) + 5(t), (1.197) 

where the normally ordered two-time correlation function for x is defined as 

{: xit + T)xit) :) = Ti- {{c + c^)e^^[cpit) + p{t)c^]} . (1.198) 

That is because the photocurrent at time t changes the systeni through the 
feedback as weü as through the conditioning. Taking this into account, it is 
not difhcult to show that the correct expression is 

E[/h°-(í + r)/h°-(í)] = í?Tr {(c + c^)e^^[{c - íF/r;)p(í) 

+ p(í)(ct +^F^)]} +(5(t). (1.199) 

Note that the feedback is present in both the term in square brackets and 
the evolution by £ for the time r. The former presence means that the in- 
loop photocurrent may have a sub-shot-noise spectrum, even if the system 
dynamics are classical. That is to say, for an optical system it may be possible 
to use an a semiclassical analysis (with classical light fields and detector shot 
noise) which nevertheless correctly predicts a sub-shot noise photocurrent 
spectrum. This is not surprising given the analogous result for free squeezing 
in Sec. 111331 

1.9 In-loop Squeezing Revisited 
1.9.1 In-Loop Squeezing 

In Sec. 11.31 it was shown that a sub-shot-noise in-loop photocurrent is not 
evidence for squeezing of the in-loop light in the usual sense, not least because 
the usual two-time commutation relations do not apply to an in-loop field. 
Also, a linear optical element (a beam splitter) fails to extract any squeezing 
from the loop, and in fact produces an above-shot-noise output. Nevertheless 
it was also shown that an in-loop nonlinear optical element (a QND intensity 
meter) agrees with the intensity statistics seen by the in-loop detector. This 
means that it is an open question as to whether the effect of the in-loop 
"squashed" light on other nonlinear optical elements is more akin to 
that of squeezed light or light with classical fluctuations. 

The simplest nonlinear optical element is a two-level atom. Shortly after 
the first observation of squeezing, Gardiner jïïü] made a seminal prediction 
regarding its effects on such an atom |60| . namely that immersing an atom in 
broad-band squeezed light would break the equality between the transverse 
decay rates for the two quadratures of the atòmic dipole. In particular, one 
decay rate could be made arbitrarily small, producing an arbitrarily narrow 
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line in the power spectrum of the atom's fluorescence. This was seen, as the 
title of inn] proclaims, as a "direct efFect of squeezing". 

In this section I present work published recently which considers 
the question of whether this atòmic hne-narrowing is characteristic only of 
squeezing in the conventional sense ('real squeezing'), or whether it can be 
produced by hght which gives rise to a below shot-noise photocurrent by 
virtue of being part of a feedback loop ('in-loop squeezing'). The answer is 
that in-loop squeezing can do the job. In fact, the inhibition of the decay of 
one quadrature depends on the amount of squeezing and the quality of mode- 
matching to the àtom in exactly the same way for in-loop squeezing as for 
free-field squeezing. Thus in-loop squeezing appears likely to be an important 
tool for future experimental investigation of the effect of low-noise light on 
àtoms, as it is usually easier to generate than free squeezing. 

Consider the apparatus shown in Fig. 11.31 but for the moment without 
the fluorescent atom. The Mach-Zehnder interferometer on the left hand side 
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Fig. 1.3. Diagram of the experimental configuration discussed in the text. AU beam 
splitters are 50:50. The atom is represented by the small ellipse at the focus of 
&in(í). The difïerence /hom(í) between the photocurrents at detectors Dl and D2 is 
amplified and split. The two signals (with opposite sign) are fed back to the two 
electro-optic modulators (EOM) 



has two functions. First it produces a weak beam (òjn) which is given by 
6i„(í) = v{t) - (í/2)(/3e*'^ - pe-'^) w i/(í) -I- (1.200) 
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Here (3, assumed real, is the coherent amplitude of the làser, and ±0, as- 
sumed small, are the phase shifts imposed by the electro-optic modulators. 
The operator ^{t) represents vacuum fluctuations as in Sec. 11.31 The second 
function of the interferometer is to produce a local oscillator beam with mean 
amplitude 

- (l/2)(/3e*'^ + /3e-'^) = -/3[1 + 0(02)]. (1.201) 

With appropriate phase shifts assumed, this local oscillator is then used for 
making a homodyne measurement of the X = 6 + quadrature of òout 
(which, in the absence of the atom, is idcntical with òin). If the cfhciency 
of the detectors is e then the homodyne photocurrent is represented by the 
operator |[2j| 



/hom(í) = %/Í^out(í) + VT^Ut), (1-202) 

where ^dt) is a unit-norm real white noise process. 

Although S'(^t(a;) is shot-noise limited for high freqüències, it need not 
be for lower freqüències. In particular, the feedback loop shown can produce 
a spectrum below the shot-noise as shown in Sec. 11.31 The current /hom(í) is 
amplified and used to control cj). If we set 



0(í) = -f^ / h{t')hon.{t -T- t')dt', (1.203) 
2PV£ Jo 

then we have a feedback loop with a low-frequency round-loop gain of g 
which is stable under the conditions given in Sec. II. 3. 21 Solving in the Fourier 
domain. 



-^out('^) 



Ç,(^)+g/ï(a;)e— ^\/^ée(t^) 



1 - gh[L·ü)e-^'^^ 



(1.204) 



Thus (which here equals Xout) has a spectrum 

S^du)^[l+g^\h{u^)\^{e-'-l)]l\l~gh{üü)e-^-^\\ (1.205) 

At a frequency íD much less than the feedback bandwidth ^ r~^, and much 
less than the reciprocal of the delay time T"^, h{üj)e^^'^'^ ~ 1 and the mini- 
mum noise is 

5ií(w)mi„ = l-£, for.g = -£/(l-e), (1.206) 



which is clearly below the Standard quantum limit. These results reproduces 
those of H1.83|l . where here e is playing the same role as 772. 
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1.9.2 An In-Loop Atom 

Returning to Fig. 11.31 we now include thc two-level atom, which is assumed 
to be resonant to the làser. It couples strongly only to modes oí the radiation 
field having the appropriate dipole spatial distribution 1 . However, by fo- 
cusing a beam as shown in Fig. 11.31 it is possible to mode-match a significant 
proportion, say 77, of 6in into the atom's input. Recent numerical calculations 
indicate that practical schemes for focusing light in frec space have a limit 
on r] of order 0.1 |62j. This suggests that in practice, a more efScient way to 
increase the efïective 77 would be to couple the light into a microcavity, as 
in 63 . However, it is conceptually simpler to consider the free-space set up 
in Fig. O 

The Hamiltonian of the atom in thc intcraction picture at time t is 

H{t) = -liVvhnit) + v/T^//(í)]at(í) + H.c. (1.207) 

Here cr = |.g) (e| is thc atòmic lowering operator and 1 have set the longitudinal 
atòmic decay ratc to unity. The operator fi{t) represents an independent vac- 
uum input. Under this coupling, the output field is found from the techniques 
of Sec. O to be 

bont{t) = bin{t) + ^ait). (1.208) 

Although it would be possible to give a description of the entire feedback 
loop in terms of atòmic and radiation field operators j51| . it is simpler to 
use the quantum trajectory theory of homodyne measurement as outlined 
in Sec. 11.61 In this theory, only the atom is treated as a quantum mechan- 
ical System with state matrix p{t)] the rest of the apparatus is considered 
as a complicated measurement and feedback device for the atom. The pho- 
tocurrent /''°'°(í) is therefore a classical quantity. From (|1.153|l it is given 

by 

I^°"\t) = ï^°'^it) + e''°'"(í), (1.209) 

where ^^°'°(í) is local-oscillator shot noise, which in this case is the only 
source of noise in the whole system. From (ll.202|) and (|1.2U8|I . the expected 
value ï^°™{t) conditioned upon the prior photocurrent record is 

I!^°^{t) = ^Tiip,{t)a.,] + V^2/3Mt)- (1-210) 

Here çic (í) is not set to its average value of zero because it is determined by 
the prior classical photocurrent via H1.203|l . 

From the theory of Sec. 11.61 the atom will obey the foUowing nonlinear 
stochastic master equation 

dp, = dW[a]p, + V^dVF^°'"(í)H[a]pe - ídí[iífb,pc]. (1.211) 

The final Hamiltonian is due to the feedback. It is identical to the term due 
to feedback in the fundamental atòmic Hamiltonian l|1.207() . namely 



(1.212) 
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So far, it is not a singular quantity because of the smoothing of the photocur- 
rent in |^ 



Now consider the limit of instantaneous feedback on the atòmic time-scale, 
T, r < 1. In this hmit 

2m) = 9lc°"'it)/V^ (1-213) 
and thus we can derive 

= (1 - g)-'{m + Vv^Tt[pS)ct.]}. (1.214) 
Hence from (ll.212|l and l|l. 213(1 . 

iífb(í) = Xlay{TT[p,{t)a,]+m/Vv^}, (1-215) 

where it is to be understood that t on the right-hand side of this equation 
actually stands for í — 0+. Since — cx) < 5 < 1, the feedback parameter A is 
given by 

A = ot/(1-5) e i-V,^)- (1-216) 

Taking the Markovian hmit ahows us to derive a deterministic master 
equation for the atom. Equation H1.211|l no longer apphes, because H^{t) in 
p. 215(1 is singular. However, by comparison of 1(1.215(1 with ((1.157(1 we see 
that we can foUow the methods of Sec. ll.G.^ and Sec Il.SI to obtain 

\2 

p = V[<t]p - i\[\cFy, ap + pa^] + —V[^(jy\p. (1.217) 

TjE 



This equation, and the foUowing relation between A and the in-loop squeezing 

, 1 + g2(e-i - 1) 2A A2 , ^ 

'5in(<^) = ;/ =1 + - + Z^^ (1-218) 



are the central results of this section. In ((1.218() . we still have ü ^ r^^, but 
now also íD 3> 1. This ensures that the atòmic variables (with characteristic 
time scale of unity) do not contribute significantly to the spectrum at w, so 
that l(1.205() is still vahd. 

From the master equation 1(1.217(1 it is easy to derive the foUowing dy- 
namical equations: 

Tr[pfT,] = -7xTr[pa,] (1.219) 
Tr[p(7j,] = -^yTilpay] (1.220) 
Tr[pa.] = -7.Tr[pa,] - C (1.221) 

Only the equation for ay is unaffected by the feedback, with 7^ — 1/2. The 
new decay rate for ax is 



1 

7. = 2 



A 



2' 

2A+ — 



(1.222) 
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and the modified parameters for are 

7í = 7a + ; C = 1 + A. 
In steady state 

Tr[pssí7j:] = TrÍpssCTy] = 0, 

Tt\p,,o,\ = -1 + AV[2r/e(l + A) + A^ 



(1.223) 



(1.224) 
(1.225) 



The most interesting of these results is that negative feedback can reduce 
the decay rate of the x component of the atòmic dipole below its natural 
value of 1/2. From l|1.218|l it can be re-expressed as 



(1.226) 



This clearly shows that has two contributions: ^(1 — 77) from the vacuum 
input and ^^riS-^^iyi) from the in-loop squeezed hght. The greatest reduction 
occurs for minimum in-loop fluctuations as in (|1.206() . for which 



i (1 — r\e) , for A 



(1.227) 



The slower decay of can be directly observed in the power spectrum 
of the fluorescence of the atom into the vacuum modes. This measures the 
photon flux per unit frequency into these modes and is defined by 



P{UJ) 



1 



27r 



^/àt(-c.)a(0))s 



From H1.219|l - 11.221|) this is easily evaluated to be 



P(^) 



(l-77)(7,-C) 
8777^ 



Ix 



là + ^ 7y + ^ 



For the optimal squeezing (A = —rje) we have 



(1 - ??)?7g 
47r(2 - -qe) 



(1 - rjs)'^ + Auj"^ 



(1.228) 



(1.229) 



(1.230) 



This is plotted in Fig.Ofor 77 = 0.8 and e = 0.95. 



1.9.3 Comparison with Free Squeezing. 

To compare the above results with those produced by free squeezing we again 
assume that the mode-matching of the squeezed modes into the atom is 77, 
and that the squeezing is broad-band compared to the atom. Assuming also 
that the input light is in a minimum-uncertainty state for the X and Y 
quadratures fB., it can be characterized by a single real number L, with 

S^^iu;) = L=l/Sliiü). (1.231) 
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Fig. 1.4. Plot of the Power Spectrum P{uj) of the fluorescence into the vac- 
uum modes, for in-loop squeezing (sòlid) and free squeezing (dotted), with mode- 
matching r; = 0.8 and squeezing S{^{lü) = 0.05. The linewidth for in-loop squeezing 
is slightly broader because the contribution from ay is not broadened in this case. 
The natural-width spectrum from weak thermal driving (dashed) is scaled up for 
comparison 

This parameter is related to those of Sec. ll.T.^ bv L — 2N + 2M + 1, where 
AP = N{N + 1) and M is assumed real. This yields the master equation 



p={l~ Tjma]p + ^V[{L + l)a - (L - 1)^^]^, (1.232) 
which leads again to H1.219|l - H1.221fl . but with 

7.^^[{1-V)+VL], (1-233) 

7y^^[{l-v)+vL-'], (1.234) 

7. -7x+7a, C=l. (1.235) 



For L < 1 the decay of is again inhibited. The crucial observation to be 
made is that the dependence of on the degree of X quadrature squeezing 
of the input light is exactly the same as for in-loop squeezing, as is seen by 
comparing H1.231|l and (|1.233|l with (|1.226|) . The only difïerence between the 
two cases is that C is unafïected by the free squeezing and that is not 
increased by the in-loop squeezing. The latter is a direct consequence of the 
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fact that an in-loop field is not bound by the usual two-time uncertainty 
relations. The free squeezing fluorescence spectrum is again given by p.229|l . 
This is also plotted in Fig.Ofor rj = 0.8 and L = 0.05. As this fi gure shows, 
the spectra are certainly not identical, but the sub-natural hnewidth is much 
the same in both. 

To conclude, hne-narrowing of an atom is not a diagnòstic of free squeez- 
ing. Rather, it requires only temporal anticorrelations of one quadrature of 
the input field (for timcs much shorter than the atòmic lifetime) such as can 
be produced by a negative electro-optic feedback loop. The dependence of the 
line-narrowing on the input squeezing and the degree of mode-matching is 
the same for in-loop squeezing as for free squeezing. Because the quadrature 
operators of an in-loop field do not obey the usual two-time commutation 
relations, the reduction in noise in one quadrature does not imply an in- 
crease in noise in the other. Hence the line-narrowing of one quadrature of 
the atòmic dipole by in-loop squeezing does not entail the line-broadening 
of the other quadrature. What significance this difference has in the physics 
of more complex atòmic interactions with squeezed light |fíO| is a question 
requiring much investigation. 

In-loop squeezing is generally easier to producc than free squeezing for a 
number of reasons. First, in-loop squeezing does not require expensive and 
delicate sources such as nonlinear crystals, but rather ofï-the-shelf electrònic 
and electro-optical equipment. Second, the amount of squeezing is limited 
only by the efíiciency of the photodetection. For homodyne detection, as re- 
quired here, an efiiciency of 95% is readily obtainable and would enable 
in-loop squeezing of 95%. Third, in-loop squeezing can be produced at any 
frequency for which a coherent source is available, so experiments could be 
conducted on any atòmic transition. The one difficulty with in-loop squeezing 
is that it requires a feedback loop responsc timc much shorter than an atòmic 
lifetime, but this would not be a problem for metastable transitions. Thus as 
well as giving us a better theoretical understanding of the efïects on matter 
of light with fluctuations below the Standard quantum limit, in-loop squeez- 
ing should be a practical alternative to free squeezing in the experimental 
investigation of these effects. 

1.9.4 Other Uses of "Squashed light" 

The analysis given above is probably the most interesting application for 
squashed light proposed so far. However, it is worth pointing out that another 
conceivable application was suggested in This was in the context of a 
QND measurement of an the squashed in-loop quadrature. For a perfect 
QND device, such a measurement does not disturb the feedback loop. Now 
because the in-loop fluctuations are squashed, the readout from the QND 
device will have less noise than it would have without feedback. This is useful 
if the measurement is being made with the aim not of determining the in- 
loop intensity, but of determining the coupling constant between the in-loop 
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light and the QND meter. For a given detector bandwidth, the squashing 
of intensity fluctuations would thus enable this parameter to be estimated 
more accurately. The improvement in accuracy would be limited only by the 
eSicieney of the photodetectors used in the feedback loop, as in (11.20611 . 

Another potential apphcation for squashed hght was proposed recently 
in '31 . Here the idea is to reduce radiation pressure fluctuations on a mir- 
ror by squashing the amphtude fluctuations. The feedback involves Standard 
photodetection and an amphtude modulator. The point of reducing the ra- 
diation pressure fluctuations is again to aUow a more accurate estimation of 
other parameters, namely the spectrum of the thermal noise in the mirror. 
In this case it turns out that in the optimal regime the fluctuations can at 
most be reduced by a factor of one half. 

1.10 Conclusion 

What then, can we say in conclusion about squeezing and feedback? The 
first and foremost fact is that, in the absence of any nonlinear optical ele- 
ments, feedback cannot produce free squeezing. A nonlinear optical element 
is any element whose effect cannot be modeled by a displacement, rotation, 
or damping of the amphtude of the light. By free squeezing I mean squeezing 
whose existence can be verified by detection in a conventional (demolition) 
detector which is not part of the feedback loop. This includes both continuum 
squeezing and intracavity squeezing (which can be measured by dumping the 
light out of the cavity). 

If there are nonlinear optical elements present then there are many in- 
teresting things one can achieve with feedback involving squeezing. First, if 
those elements are used to generate a squeezed beam which is split, then 
feedback (or indeed feedforward) onto an intensity modulator can effectively 
transfer amphtude squeezing electronically from the beam which is detected 
to the free beam. This was discussed in Sec. 11.3.71 Second, nonlinear optical 
elements can create quantum correlations between two beams. Then, even 
if the beams separately are not squeezed, feedback (or, equally effectively, 
feedforward) of the light from one beam can make the other beam squeezed. 
If they are squeezed, then feedback or feedforward can enhance this squeez- 
ing. A number of scenarios were discussed in Sec. 11.51 Third, if the nonlinear 
element creates squeezing within a cavity then controlling the driving of the 
cavity by a photocurrent derived from a homodyne measurement of the out- 
put can enhance the squeezing, as discussed in Sec. II. 71 This could be seen as 
a generalization of the first case discussed above, where the cavity mirror is 
acting as a beam splitter and the intracavity light is being continually split. 

A quite separate issue is the nature of the light within a feedback loop. At 
first sight it appears easy to make this squeezed, as negative feedback (with- 
out the aid of any optical nonlinearity) can produce an in-loop photocurrent 
with arbitrarily low noise, regardless of the efficiency of the detector. How- 
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ever, the fact that this is the current from an in-loop detector means that 
a sub-shot-noise spectrum does not have the usual significance. First, there 
is a perfectly good semiclassical explanation for this phenomenon in terms 
of coherent field states. Second, the usual two-timc commutation relations 
for a free field do not hold for an in-loop ficld. This means that there is no 
"squeezing" of the uncertainty fi'om the amplitude to the phase quadrature, 
rather just an apparent "squashing" of the amplitude uncertainty with no 
effect on the phase uncertainty. Lastly, from a practical point of view, the 
squeezed light cannot be removed from the loop using a linear optical device. 
A beam splitter inserted in the in-loop beam produces a free (out-of-loop) 
beam with noise level above rather than below the shot noise. 

Despite the differences between in-loop squashing and free squeezing, it 
turns out that there are similarities. An in-loop quantum non-demolition 
amplitude detector will respond to amplitude squashing exactly as to am- 
plitude squeezing. For a feedback loop constructed using a perfectly efficient 
(demolition) detection, a perfect in-loop non-demolition detector will have a 
read-out identical to that of the in-loop demolition detector, and, for arbi- 
trarily large negative feedback, this will have an arbitrarily low noise level. 
For non-unit efficiency of the demolition detector in the feedback loop, the 
two detectors will not agree and the maximum degree of squashing observed 
by the non-demolition detector will be limited by the demolition detector 
efficiency. 

The degree of squashing seen by a perfect QND detector turns out to 
be a legitimate measure of the degree of noise reduction in the light, and 
this squashing can have effects on nonlinear optical devices very similar to 
those produced by squeezing. In particular, shining broad-band quadrature 
squashed light on a resonant two-level atom will cause the decay of the in- 
phase quadrature of the atòmic dipole to be suppressed. This manifest itself 
as a line-narrowing in the power spectrum of the atom's fluorescence, an effect 
which was originally thought to be characteristic of squeezing. Moreover, the 
reduction in the linewidth of the in-phase quadrature depends on the degree 
of squashing (as measured by a QND device) in precisely the same way as it 
does on the degree of squeezing in the case of freely propagating light. This 
has important experimental implications as there are are a number of factors 
which make squashing easier to achieve than squeezing. 

Very recently, there have been more intriguing developments in the theory 
of squashed light and its application to quantum spectroscopy "ïïïï". First, I 
showed that it is possible to squash light simultaneously in both quadratures, 
and also to squeeze light in one quadrature and squash it the other. In the 
limit of perfect squeezing, and perfect squashing (which requires unit effi- 
ciency detectors) the in-loop fluctuations can be banished from both quadra- 
tures! An atom coupled only to this in-loop light would (to a first approxi- 
mation) have its spontaneous decay completely suppressed. This surprising 
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prediction is one more example of the continuing fruitful investigation of the 
relationship between feedback and squeezing. 
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